Deep ROSAT-HRI observations of the NGC 1399/NGC 1404 
region: morphology and structure of the X-ray halo 
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ABSTRACT 

^ We present the analysis of a deep (167 ks) ROSAT HRI observation of the cD galaxy NGC 1399 

in the Fornax cluster, comparing it with previous work on this galaxy and with recent Chandra 
data. We find, in agreement with previous observations, an extended and asymmetric gaseous 
halo with a luminosity (in the 0.1-2.4 keV energy band) oiLx = (5.50±0.04) x 10**^ erg s~^ within 
46 kpc (assuming a distance of D=19 Mpc). Using both HRI and, at larger radii, archival PSPC 
^> ' data, we find that the radial behavior of the X-ray surface brightness profile is not consistent with 

, a simple Beta model and suggests instead three distinct components. We use a multi-component 

' bidimensional model to study in detail these three components that we identify respectively with 

the cooling flow region, the galactic and the cluster halo. From these data we derive a binding 
' mass distribution in agreement with that suggested by optical dynamical indicators, with an 

inner core dominated by luminous matter and an extended dark halo differently distributed on 
galactic and cluster scales. The HRI data and a preliminary analysis of Chandra public data, 
\ allow us to detect significant density fluctuations in the halo. We discuss possible non-equilibrium 

fS i' scenarios to explain the hot halo structure, including tidal interactions with neighboring galaxies, 

ram stripping from the intra-cluster medium and merging events. 
\^ • In the innermost region of NGC 1399, the comparison between the X-ray and radio emission 

0^ \ suggests that the radio emitting plasma is displacing and producing shocks in the hot X-ray 

' emitting gas. We do not detect the nuclear source in X-rays and we pose an upper limit of 

~ 4 X 10^^ erg s~^ (0.1-2.4 keV) to its X-ray luminosity. 

We found that the NGC 1404 halo is well represented by a single symmetric Beta model 
^ , and follows the stellar light profile within the inner 8 kpc. The mass distribution is similar to 

■ the 'central' component of the NGC 1399 halo. At larger radii ram pressure stripping from the 

intra-cluster medium produces strong asymmetries in the gas distribution. 

Finally we discuss the properties of the point source population finding evidence of correlation 
between the source excess and NGC 1399. 

Subject headings: X-rays: individual (NGC 1399, NGC 1404) — galaxies; halos — galaxies: clusters: indi- 
vidual (Fornax) — galaxies: jets — radio continuum: galaxies 

1. INTRODUCTION 
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extensively studied in a wide range of wavelengths, 
ranging form radio to X-rays. The optical radial 
profile, first studied by Schombert (1986) and later 



1 



in more detail by Killeen & Bicknell (1988), re- 
veals a large halo extending up to 250 kpc. The 
galaxy is surrounded by a large number of glob- 
ular clusters, 10 times in excess with respect to 
those of its nearer companion NGC 1404 and of 
the Fornax galaxy NGC 1380 (Kissler-Patig et al. 
1999). Dynamical studies of the stellar popula- 
tion of NGC 1399 indicate that the central regions 
of the galaxy are characterized by a high velocity 
dispersion, decreasing with radius, and slow ro- 
tation, as expected in giant ellipticals (Graham 
et al. 1998). At larger radii the velocity disper- 
sion starts increasing again, as shown by studies 
of the planetary nebulae (Arnaboldi et al. 1994) 
and globular clusters dynamical properties (Grill- 
mair et al. 1994). This suggests a different dynam- 
ical structure, than that of the inner stellar body, 
for the galaxy envelope and the globular cluster 
population. 

NGC 1399 hosts a weak nuclear radio source 
(Killeen, Bicknell & Ekers 1988) with radio lumi- 
nosity of ~ 10^^ ergs s^^ between 10^ and 10^^° 
Hz. The radio source has two jets ending in dif- 
fuse lobes, confined in projection within the opti- 
cal galaxy. Killeen, Bicknell & Ekers (1988) sug- 
gest that the radio source is intrinsically small and 
confined by the thermal pressure of the hot ISM. 

X-ray data have shown the presence of an ex- 
tended hot gaseous halo surrounding NGC 1399. 
Observations made with the Einstein IPC by Kim, 
Fabbiano & Trinchieri (1992) constrained the gas 
temperature to be kT>l.l keV. The Einstein IPC 
data were used also by Killeen & Bicknell (1988) to 
derive the mass distribution of the galaxy. They 
found that both models with and without dark 
matter were compatible with those data, depend- 
ing on the assumed temperature profile. White 

(1992) suggested the presence of a cooling flow in 
the center of the galaxy, depositing 0.8±0.6 Mq 
y~^. Ginga observations (Ikebo et al. 1992) led to 
the detection of extended emission out to a galac- 
toccntric radius of ~250 kpc. Serlemitsos et al. 

(1993) constrained the X-ray temperature to be 
1.0<kT<1.5 keV with the Broad Band X-ray Tele- 
scope. Rangarajan et al. (1995, hereafter RFFJ), 
using the Rosat PSPC, have studied in detail the 
temperature profile of the hot inter-stellar medium 
(ISM) out to 220 kpc finding an isothermal profile 
(0.9<kT<l.l) from 7 to 220 kpc, and a central 
cooling flow (0.6<kT<0.85) of at least 2 Mq y-^. 



Ikebe et al. (1996) were able to identify with 
ASCA the presence of different components in 
the X-ray halo, associated respectively with the 
galaxy and the Fornax cluster potential. Jones 
et al. (1997) used the better resolution of the 
Rosat PSPC for studying in detail the cluster X- 
ray structure and flnding a total binding mass of 
(4.3-8.1) X 10^2 Mq within 100 kpc, with a mass- 
to-light ratio increasing from 33 Mq /Lq at 14 kpc 
to 70 M0/L0 at 85 kpc. They also found a metal 
abundance of 0.6 solar. 

Buote (1999) analized ASCA data finding that 
either a two temperature spectral model or a 
Cooling-Flow model with Fe solar abundances, are 
required to fit the thermal X-ray emission, and ad- 
ditional absorption is needed in the galaxy center. 
Matsushita, Ohashi & Makishima (2000) found 
similar near solar metallicities but argued that the 
metallicity inferred by ASCA spectral fits is de- 
pendent on the assumed atomic physics model and 
is not solved by muti-components spectral models. 

More recently Sulkanen & Bregman (2001) used 
ROSAT data to pose an upper limit on the nuclear 
source brightness. The Chandra data were used 
by Loewenstein et al. (2001) to further constrain 
the nuclear source and by Angelini, Loewenstein & 
Mushotzky (2001) to study the point source pop- 
ulation hosted by globular clusters. 

In this paper we present the results of a deep 
image of the NGC 1399-NGC 1404 field obtained 
from data collected between 1993 and 1996 with 
the ROSAT High Resolution Imager (HRI, for a 
description see David et al. 1996). We take advan- 
tage of the 5" resolution of the HRI to study in 
detail the structure of the galactic halo and relate 
the results to those obtained at larger scales with 
poorer resolution instruments. We study the in- 
teractions between the nuclear radio source and 
the galactic halo and discuss the properties of the 
discrete sources population. A preliminary analy- 
sis of Chandra data supports the Rosat results. 
Throughout this paper we adopt Hq = 75 km s~^ 
Mpc~^ and a distance of 19 Mpc (l'=5.5 kpc). 

2. OBSERVATIONS AND DATA ANAL- 
YSIS 

The NGC 1399 field, including NGC 1404, 

was observed at three separate times with the 
ROSAT HRI: in February 1993, between January 
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Table 1: Rosat HRI observations of NGC 1399. 



name 


Field center 


sequence id. 


Exp. time 


obs. date 


P.I. 




R.A. 


Dec 




(sec) 






NGC 1399 


03'^38™31= 


-35°27'00" 


RH600256n00 


7265 


1993 Feb 17 


D.-W. Kim 




)) 


)! 


RH600831n00 


72720 


1996 Jan 04-1996 Feb 23 


G. Fabbiano 




55 


5) 


RH600831a01 


87582 


1996 Jul 07-1996 Aug 26 


G. Fabbiano 



and February 1996 and between July and Au- 
gust of the same year. The total exposure time 
is 167.6 Ks (Table 1). The data were processed 
with the SASS7.8 and SASS7.9 versions of the 
ROSAT standard analysis software (SASS). For 
our data analysis, we used the IRAF/XRAY and 
CIAO packages developed at the Smithsonian As- 
trophysical Observatory and at the Chandra X-ray 
Center (CXC), and other specific software as men- 
tioned in the text. 

2.1. Aspect Correction 

The HRI data processed with the SASS versions 
prior to SASS7_B of March 1999 (as it happens for 
our data) suffer from an error in the aspect time 
behavior (Harris 1999). This translates into an er- 
ror on the position of the incoming photons and 
thus it results in a degradation of the Point Re- 
sponse Function (PRF). Because high resolution 
analysis is the primary objective of this work we 
run the correction routine ASPTIME (F. Primini 
2000, private communication) on the data, to im- 
prove the aspect solution. Visual inspection of the 
brightest pointlike sources in our field (Figure 1) 
demonstrates the improvement in the image qual- 
ity. 

Additional problems that may preclude the at- 
tainment of the HRI potential resolution are the 
imperfect correction of the spacecraft wobble and 
the wrong tracking of reference stars due to the 
variable pixel sensitivity across the detector (Har- 
ris et al. 1998). We followed the procedure sug- 
gested by Harris and collaborators, of dividing 
the observations in time bins (OBI) and realigning 
these segments by using the centroids of a refer- 
ence pointlike source in each OBI. 

Our data can be divided in 4 OBIs in RH600256n00, 
43 in RH600831n00 and 33 in RH600831a01. We 
checked each OBI individually using the bright- 
est point source in the field (No. 8 in Table 3). We 
found that in the first observation (RH600256n00) 




Fig. 1. — Comparison between the aspect of the 
brightest pointlikc source in the NGC 1399/1404 
field (No.8 in Table 3) before (left) and after 
(right) correction. The pixel size is 4". Contour 
levels are spaced by 5.9 x 10^^ cnts arcmin^^ s^^ 
with the lower one at 6 x 10~^ cnts arcmin"^ s~^. 

the OBI are well aligned and no correction is 
needed. The second observation (RH600831n00) 
has 7 OBIs out of 43 that are 7"-10" displaced; 
however, the overall PRF is 6" FWHM so that just 
a slight improvement can be obtained dcwobbling 
the image. The third observation (RH600831a01) 
has a 8" PRF but the signal to noise ratio of the 
individual OBIs is worse than in the previous ob- 
servations so that the dewobbling procedure can- 
not be applied. Given these results, we decided 
not to apply the dewobbling correction, and to 
take as good a resulting PRF, in the composite 
image, of ~7" FWHM. 

2.2. Composite Observation and Exposure 
Corrections 

The aspect-corrected observations were co- 
added to obtain a "composite observation". To 
make sure that pointing uncertainties did not de- 
grade the image we used the centroids of three 
bright pointlike sources in the field (No. 8, 19 an 
24 in Table 3) to align the images. The applied 
corrections were all within 2 arcsec. The resulting 
composite image is shown in Figure 2. This image 



3 



I I I I I 




■ I I I I I I I ■ 

3:40:00 39:00 38:00 00 

Right Ascension (J2000) 

Fig. 2.— The 'raw' NGC 1399/NGC 1404 HRI field. The composite image is here displayed after it was 
rebinned to 5 arcsec/pixel. Even without exposure correction the extended emission surrounding NGC 1399 
is clearly visible. The elongated features in the lower left corner are due to the presence of "hot spots" on 
the detector. 



was then corrected for vignetting and variations of 
exposure time and quantum efficiency across the 
detector by producing an "exposure map" for each 
observation with the software developed by S.L. 
Snowden (Snowden et al. 1994, hereafter SMB). 

2.3. Brightness Distribution and X-ray/Optical 
Compeirison 

To study the large-scale brightness distribution 
of the NGC 1399 field we rebinned the exposure 
corrected data in 5"x5" pixels. We then adap- 
tively smoothed the image with the CXC CIAO 
csmooth algorithm which convolves the data with 
a gaussian of variable width (depending on the 
local signal to noise range of the image) so to en- 
hance both small and large scale structures. 

The resulting image (Figure 3a) shows a com- 
plex X-ray morphology. The center of the image is 



occupied by the extended halo of NGC 1399. The 

galaxy possess a central emission peak and an ex- 
ternal extended and asymmetric halo. This halo 
is not azymuthally symmetric with respect to the 
X-ray peak, but it extends more on the SW side. 
The X-ray surface brightness distribution of the 
halo appears filamentary, with elongated struc- 
tures and voids. As it can be seen from a compar- 
ison with the optical Digitized Sky Survey (DSS) 
image (Figure 3b), while the X-ray emission peak 
is centered on the optical galaxy, the X-ray halo 
extends radially much further than the optical dis- 
tribution. Moreover most of the features seen in 
the X-ray image have no direct optical counter- 
part. 

The X-ray emission of NGC 1404 instead is al- 
most symmetric and largely consistent with the 
optical distribution, with the exception of the SE 
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Fig. 3. — (a) 5 X 5 arcsec/pixel adaptively smoothed image of the central part of the NGC 1399/1404 HRI 
field. Colors from black to yellow represent logarithmic X-ray intensities from 5.6 x 10^'^ to 3.58 x 10^^ cnts 
arcmin"^ s~^.(b) X-ray brightness contours overlaid on the 1 arcsec/pixel DSS image (logarithmic grayscale). 
Contours are spaced by a factor 1.1 with the lowest one at 6.1 x 10~^ cnts arcmin"^ s~^. The X-ray emission 
peak is centered on the optical galaxy for both NGC 1399 and NGC 1404. In the case of NGC 1404 the 
X-ray isophotes are consistent with the optical distribution while in NGC 1399 the X-ray emission extends 
further out than the optical one. 



tail reported by Jones et al. (1997) in the PSPC 
data, that may be the signature of ram pressure 
stripping of the NGC 1404 corona infalling toward 
the denser NGC 1399 cluster halo. 

A number of pointlike source in the field is also 
evident in Figure 3a. The analysis of these sources 
is discussed in § 2.9. 

We examined the inner halo of NGC 1399 in 
greater detail using a 1" /pixel resolution. The 
adaptively smoothed image (Figure 4) reveals an 
elongation of the inner halo structure in the N-S 
direction plus a large arc protruding to the West 
side of the halo. Several voids are present in the X- 
ray distribution, the largest being 1.5 arcmin NW 
of the emission peak. The central peak appears 
circular with at most a slight N-S elongation. 

These features are all above 3(7 significance, 
since 3cr was the minimum significance level 
requested for an intensity fluctuations to be 
smoothed on a given scale. Some features, as 
the SW clump centered at RA, Dec=3''38'°15%- 



35°29'00" or the one at 3'^39'"10",-35°26'30" , are 
also detected with a wavelets detection algorithm 
(see § 2.9). We can rule out that features on scales 
of 1 arcmin are of statistical nature, because their 
significance can be as high as 6cr (see § 2.5). These 
brightness fluctuations must be due to the local 
physical conditions of the hot gaseous halo. 

2.4. Radial Brightness Profiles 

24.1. NGC 1399 

As a first step in the quantitative study of 
the X-ray emission, we created a radial profile 
from the HRI data, assuming circular symmetry. 
Count rates from the exposure-corrected compos- 
ite image, were extracted in circular annuli cen- 
tered on the X-ray centroid RA, Dec=3''38'"28".9, 
-35°27"02'.l. We took care to remove the con- 
tribution of all the detected point-like sources by 
excluding circles within the 3t7 radius measured 
by the wavelets algorithm (§ 2.9), from the source 
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Fig. 4. — Adaptively smoothed image of the inner 
NGC 1399 halo. Colors from black to yellow repre- 
sent logarithmic X-ray intensities from 6.4 x 10^'^ 
to 7.7 X 10^"'^ cnts arcmin^^ s^"'^. 

centroid. We also excluded a 150" circle centered 
on NGC 1404 (see § 2.4.2). 

The X-ray brightness profile is shown as a con- 
tinuous line in Figure 5. The emission extends out 
to ~ 500" (46 kpc at the assumed distance of 19 
Mpc). The radial profile flattens out past 500" 
(dash-dotted line in Figure 5) suggesting that we 
have reached the field background level and that 
any residual halo emission is below our sensitiv- 
ity limit. Nevertheless we know from previous 
investigations of the Fornax cluster (Jones et al. 
1997; Ikebe et al. 1992; Killeen & BickneU 1988) 
that extended X-ray emission, bright enough to be 
detected in our data, is present at galactocentric 
radii > 500". We thus expected to see a smooth 
gradient in the X-ray profile continuing past 500" . 
We checked to see if a very high background was 
present in some OBIs or either in some of the HRI 
spectral channels but we found none that could 
explain our lack of sensitivity. A possible expla- 
nation is that we didn't reach the background level 
around 500" and the actual background is much 
lower but, due to uncertainties in the exposure cor- 
rection near the edge of HRI field of view, we are 
not able to see the expected decline in the X-ray 
emission. 



1.000 




r (arcsec) 

Fig. 5.— HRI radial profile of NGC 1399 X-ray 
surface brightness (continuous line) after exclud- 
ing aU point sources (§ 2.9) and NGC 1404. The 
dot-dashed line represents the flattening level of 
the HRI profile measured in the 500" -850" annu- 
lus. Short and long dashed lines represent respec- 
tively the SMB background level before and after 
rescaling to match PSPC counts. The rescaled 
background was used to derive the background- 
subtracted counts (open circles) and the best-fit 
Beta model within 50" shown in Figure 8 (dot- 
ted line). The HRI PRF range is represented by 
the shaded region. Radial bins are 2" wide up to 
r=50" , and 5" wide at larger radii. 

We used the SMB software to calculate a back- 
ground map for this field. The radial profile de- 
rived from this background map is shown as a 
small dashed line in Figure 5. The SMB model 
however may underestimate the "true" HRI back- 
ground somewhat because it models only the 
charged particle contribution. We therefore cor- 
rected this value using the 52 ks Rosat PSPC ob- 
servation RP600043N00 centered on NGC 1399, 
taken on the 1991 August 15 (Figure 6). Even if 
the PSPC observation is much shorter than our 
HRI image, the higher sensitivity and lower back- 
ground of this instrument allows a better study of 
the large scale X-ray emission. 

The exposure corrected PSPC profile is shown 
in Figure 7. The background level (dashed line) 
was estimated from the 2400" -3000" annulus. We 
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Fig. 6. — Exposure corrected PSPC image of the 
NGC 1399/1404 field. The two brightest central 
emission peaks represent the two galaxies (NGC 
1399 and NGC 1404) visible also in the HRI field. 
X-ray contours, spaced by a factor of 1.3 with the 
lowest level at 3.5x 10~^ cuts arcmin"^ s~^, clearly 
show the presence of an extended halo much larger 
than the HRI FOV (central 40 x 40 arcmin). Ra- 
dial ribs are due to the PSPC support structure. 

compared our radial profile to the one derived by 
Jones et al. (1997) and RFFJ, who made use of 
the same data, finding consistent results. We then 
rescaled the SMB background to find the best 
agreement between HRI and PSPC counts in the 
50-400 arcsec region, after rebinning the HRI data 
to match the wider PSPC PRF. The "true" HRI 
background level found in this way (long-dashed 
line in Figure 5) lies in between the flattening level 
of the HRI profile and the SMB background and 
was used to derive the HRI background-subtracted 
profile shown as empty circles. This way of find- 
ing the HRI background assumes that the PSPC 
brightness profile is not affected by systematic er- 
rors on the background determination and can 
thus be used to rescale the HRI counts. However, 
even if the PSPC background suffered from resid- 
ual errors (e.g. due to uncertainties in the expo- 
sure correction at large radii), they would not sig- 
nificantly aff'ect the central region profile - used 
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Fig. 7. — Exposure corrected PSPC radial pro- 
file of NGC 1399 (continuous line). The dashed 
line represents the background level measured in 
the 2400-3000 arcsec region, used to derive the 
background-subtracted counts (open circles). The 
X-ray halo clearly extends at much larger radii 
than visible in the HRI FOV but the central re- 
gion is not well resolved due to the larger PRF. 

to rescale the SMB HRI background - because 
there the X-ray brightness is one order of mag- 
nitude larger than the PSPC background level. 

The HRI radial surface brightness profile has 
a complex structure. In the 'central' region, i.e. 
within 50" from the emission peak, the X-ray emis- 
sion is well represented by a simple Beta model: 

-3/3+0.5 



S oc [1 + {r/roY 



(1) 



with best fit parameters tq = 3.93 ± 0.16 arcsec, 
(3 = 0.506 ± 0.003 and = 14.5 for 22 degrees 
of freedom (Figure 8). To determine the best 
fit parameters the model was convolved with the 
HRI on-axis PRF (David et al. 1996). At radn 
larger than 1 arcmin the surface brightness pro- 
file shows a significant excess over this model (see 
Figure 5), indicating the presence of an additional 
'galactic' component. Restricting the fit to r < 
40" to check if this more extended component 
contributes around 50", gives consistent results 
with ro = 3.88tg:i^ arcsec and f3 = O.bMtoml 
(x^ = 11 for 17 degrees of freedom). Thus, 
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Fig. 9. — Background subtracted radial profile of NGC 1399 in the NE (a) and the SW (b) sectors. Open 
circles represent HRI data (2" annuli for r < 50", 5" for r > 50"), while filled triangles represent PSPC data 
(15" annuli for r < 300", 30" for r > 300"). PSPC counts arc reduced by a factor - 0.3 to match the HRI 
profile; the inner 30" are not shown because of the lower resolution. The dashed line represents the optical 
profile (convolved with the HRI PRF) from Killeen & Bicknell (1988). 



for 3" < r < 50" the X-ray emission falls as 
j,-2.04±o.02 pg^jj seen by the comparison 

with the HRI PRF (shaded region in Figure 5), 
whose uncertainty is due to residual errors in the 
aspect solution (David et al. 1996), the central 
component is extended and cannot be due to the 
presence; of a nuclear point source, that in fact is 
not detected in our data (sec § 2.8). 

The PSPC profile in Figure 7 shows that X-ray 
emission extends to much larger radii than visi- 
ble in the HRI data and suggests the presence of 
a third component dominant at r > 400" . Large 
azimuthal differences are seen in the X-ray pro- 
file at radii larger ^ 1 arcmin. We derived com- 
posite HRI-PSPC X-ray profiles in two sectors di- 
viding the images along the line connecting the 
X-ray centroids of NGC 1399 and NGC 1404: a 
NE sector (P.A. 331.5°-151.5°) and a SW sector 
(P. A. 151.5°-331.5°). While the emission within 
1 arcmin is azimuthally symmetric, the galactic 
component (r > 1') is much more extended in 
the SW direction and presents a steep decline in 
the NE sector, as suggested by Figures 3a and 4. 



The even more extended component visible in the 
PSPC data follows the opposite behavior, being 3 
times brighter around 1000" on the NE side than 
on the SW one. While these asymmetries have 
been reported by Jones et al. (1997), their PSPC 
data alone do not have enough resolution to clearly 
separate the central component from the external 
halo. Jones and collaborators used a single power- 
law model to fit the whole distribution, thus ob- 
taining a value of -^0.35 for the power law slope. 

2.4.2. NGC 1404 

We derived the brightness profile of NGC 1404 
extracting count rates in 2" annuli centered on 
the galaxy centroid. The background determina- 
tion for this galaxy is even more difficult than for 
NGC 1399, due to the fact that it is embedded 
in the halo of the dominant galaxy. Thus we had 
no other choice than to estimate the background 
level from the region between 200" and 300" from 
the galaxy (see Figure 10a), obtaining 5.60 x 10~^ 
counts arcmin"^ sec~^. The measured count rates 
are shown in Figure 10a as a continuous line while 
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(a) (b) 

Fig. 10. — (a) Radial profile of NGC 1404 extracted in 2" annuli (continuous line). Circles represent the 
background subtracted profile using the count rate measured in the 200" -300" interval (dashed line). The 
optical profile from Forbes et al. (1998) is shown as a dashed line. The HRI on-axis PRF range is represented 
by the shaded region, (b) Radial profiles extracted in 5" annuli in the North-West (P.A. 270°-360°, thin 
dashed line) and the South-East (P.A. 90°-180°, continuous line) sectors. The optical profiles and background 
levels are represented respectively by the thick dashed line and the dotted line. 



the background subtracted profile is represented 

by open circles. The NGC 1404 emission extends 
out to 110" (~10 kpc) from the galaxy center and 
shows no sign of a multi component structure as 
the one detected in NGC 1399. The surface bright- 
ness is circularly symmetric (cf. Figure 3a) and 
declines with a power law profile out to 90" . Past 
this radius the effect of ram pressure stripping be- 
comes evident as a sharp cutoff in the radial pro- 
file. In Figure 10b we can see that in the North- 
West sector (thin dashed line) the galactic halo 
suddenly disappears in the X-ray background for 
r > 90" (r > 8 kpc) while in the South-East one 
(continuous line) the emission from the tail, seen 
in the adaptively smoothed images (Figure 3), ex- 
tends up to ~ 150" (14 kpc). 

Wc fitted the observed brightness distribution, 
within a 90" radius, with a Beta model. In this 
case the on-axis PRF we used for NGC 1399 is no 
longer valid because NGC 1404 is located 10 ar- 
cmin South- West of NGC 1399. Although an an- 
alytical correction of the on-axis PRF is possible 
(David et al. 1996), it doesn't take into account the 



PRF azimuthal asymmetry. We thus convolved 

the Beta model with the on-axis PRF and consid- 
ered the resulting core radius as an upper limit to 
the correct value, obtaining ro = 6.2ltQ;2^ arcsec 
and P = 0.513to with = 51.5 for 42 d.o.f. 
The best fit model is shown in Figure 11. 

2.5. Bidimensional Halo Models 

Because of the complex halo structure it is not 
possible to obtain a satisfactory fit of the en- 
tire surface brightness distribution of NGC 1399 
with a single circularly symmetric model. The 
Sherpa modeling and fitting application of the 
CIAO package allows us to adopt a more realistic 
approach by performing a two dimensional fit of 
the X-ray images, using a complex model obtained 
by adding together three bidimensional Beta pro- 
files of the form: 

(l + (^/^0)2)3/3-0.5 (2) 
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Fig. 8. — The best-fit model and the residuals of 
NGC 1399 brightness profile within 50" . The 66%, 
90% and 99% contour levels relative to the core 
radius ro and slope /3 are shown in the inner panel. 



Fig. 11. — The best-fit model and residuals of 
NGC 1404 brightness profile within 90" . The 66%, 
90% and 99% contour levels relative to the core ra- 
dius ro and slope /3 are shown in the inner panel. 



where 

r{x,y) = "^"'^ ^ '^/^^ynew 

and 

Xnew = {x- xq) cos{6) + {y - Vq) s\n{9) (4) 

Vnew = {y- Vo) C0S{9) - {x - Xq) Sm{e). (5) 

Expression (2) represents the bidimensional exten- 
sion of the Beta model in equation (1) where we 
also allow the model ellipticity (e) and position 
angle (6) to vary. 

To model simultaneously the small and large 
scale components we used the following procedure: 
(a) we rebinned the HRI image in 15" xl5" pixels, 
in order to have enough counts per pixel to fit the 
galaxy halo and to smooth out small scale asym- 
metries; NGC 1404 and all the detected point- 
like sources (§ 2.9) were masked out before per- 
forming the fit. (b) The PSPC data were re- 
binned in 45" x45" pixels; the regions contain- 
ing the support structure where masked, together 



with all the sources detected in the HRI field. Be- 
cause the PSPC field of view (FOV) is larger than 
the HRI one, to mask all the remaining PSPC 
sources we used the source parameters measured 
by the wavelets algorithm in the context of the 
GALPIPE project (Mackie et al. 1996, and ref- 
erences therein), which has produced a list of all 
sources present in the ROSAT PSPC archive^, (c) 
We built a model composed of three bidimensional 
(3 components represented by expression (2) plus 
a constant background, (d) The first two compo- 
nents, representing the central and galactic halo 
emission, were fitted to the 15"xl5" HRI data, 
fixing the third component parameters to some 
initial values. Because at this resolution it is not 
possible to properly determine all the parameters 
of the central component, we fixed the core radius 
to the best-fit value obtained from the radial fit 
shown in Figure 5. While the central component 
shows a slight N-S elongation (e ~ 0.05, §2.3) at 

•^An online version of the GALPIPE database and docu- 
mentation can be found at the D.I. A.N. A. homepage at 
the Palermo Observatory: http://dbms.astropa.unipa.it/ 
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Table 2: Best-fit parameters for the bidimensional multi-component halo model. 



Component 


Center Position 




P 


e 


e 




V 




R.A. 


Dec 


(arcsec) 




/--I minor axis \ 
^ maior axis ' 


(rad) 




(d.o.f.) 


Central 


03'^38™28''.9 


-35°27'01" 


3.93^ 


0.54±0.02 


0.0^ 


0.0^ 


1.2" 


1612 


Galactic 


03''38™25=.0 


-35°27'40" 


2125± 825 


41± 25 


0.018± 0.018 




1.2b 


1612 


Cluster 


03''38™48=.2 


-35°23'05" 


9047± 495 


24± 3 


0.14± 0.02 


0.0 ± 0.1 


1.5 


4790 



"The Central component parameters ro, e and Q have no error because they were held fixed during the fit. 

''The value is the same for the 'Central' and 'Galactic' components because they were fitted together on the HRI image. 



this resolution our attempts to vary also ellipticity 

e and position angle B resulted in no improvements 
of the fit, so that they were fixed to zero, (e) 
The third component, corresponding to the cluster 
emission, was fitted to the 45"x45" PSPC data, 
fixing the central and galactic halo parameters to 
the values obtained in step (d). The background 
value was fixed to the flattening level measured 
in §2.4.1. (f) Steps (d) and (e) were repeated it- 
eratively until the best fit parameters converged 
within the errors. 

The best-fit model parameters are shown in Ta- 
ble 2. In Figure 12 we show the model contours 
superimposed on the 45" /pixel PSPC image (main 
panel); the nuclear region is shown in greater de- 
tail in the right-bottom panel, with contours over- 
laid on the 15" /pixel HRI image. The three com- 
ponents have different features: while the central 
one is centered on the galaxy and is highly peaked, 
the more extended galactic halo is displaced to 
the South- West thus accounting of the observed 
asymmetries shown in Figure 9. In contrast with 
the circular symmetry of the smaller components 
{.^Central ~ 0.05, ^Galactic = 0.02) the third Com- 
ponent is elongated in the East- West direction and 
its center is displaced to the North-East of the op- 
tical galaxy. These results suggest a different ori- 
gin of each component, that will be discussed in 
detail in § 3. 

The values obtained from the fit are quite 
large. This is not surprising considering the large 
spatial anisotropy of the X-ray surface bright- 
ness (see § 2.3). The situation is worse for the 
PSPC data (sec Figure 6), also because of the 
poor knowledge of the exposure corrections past 
~ 1100" from the center of the FOV (the posi- 
tion of the support ring) and the large and irreg- 
ular PRE at these radii that makes difficult a cor- 



rect source masking. In Table 2 the 'Central' and 
'Galactic' components have the same value be- 
cause they were fitted together to the HRI data. 
If we restrict the fit to the region within 1 arcmin 
and use higher resolution (2" /pixel) we obtain a 
X^ ~ 0.6, in agreement with the results obtained 
in § 2.4.1 from the profile fitting. 

The relative contribution of each component to 
the global profile is shown in Figure 13. We must 
notice that here both global and individual com- 
ponents profiles were extracted in annuli centered 
on the X-ray centroid (rather than centering each 
one on the respective component centroid) so to 
allow a comparison of the observed data with the 
model. The bottom panel shows the radial resid- 
uals respect to the best fit model. The high x^ 
values are in agreement with the values obtained 
from the bidimensional fit; in particular the resid- 
uals are systematically higher at radii larger than 
1100" due to the problems explained above. 

The residuals of the bidimensional fit to the 
HRI image are shown in Figure 14. The presence 
of structures in the gaseous halo is evident and 
confirms the results obtained with the adaptive 
smoothing technique. The highest residuals are 
coincident with the positions of point soiirc;c;s de- 
tected in the field, as can be seen from the overlaid 
contours. Considering the fiuctuations on scales of 
60" (~ 5 kpc) the statistical errors are reduced of a 
factor 4 with respect to the grayscale levels shown 
in Figure 14, and the significance of the residual 
emission is increased of the same factor. This 
means, for instance, that the deep 'hole' on the 
North- West side of the halo is approximately sig- 
nificant at the 6cr level above the fluctuations ex- 
pected from the multicomponent model, while the 
excess emission regions on the North and South 
sides of the galactic center are significant at the 
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Fig. 12. — Main panel: contours of the bidimcnsional halo model superimposed on the 40" /pixel PSPC 
image. The cluster component is clearly displaced with respect to the galactic halo; the central component 
is not visible at this resolution. Contours are spaced by a factor of 1.1 with the lowest one at 3.1x10"^ cnts 
arcmin"^ s~^. Bottom-right panel: bidimcnsional model contours overlaid on the 15" /pixel HRI image, 
showing in greater detail the galactic halo region. Contours are spaced by a factor of 1.1 with the lowest one 
at 6.0x10"^ cnts arcmin"^ s~^. 



3cr level. Even though the multi-component bidi- 
mcnsional model may not be a proper physical rep- 
resentation of the hot X-ray halo this result shows 
that the filamentary structures seen in the halo 
are not due to statistical fluctuation over a smooth 
surface brightness distribution. This result is fur- 
ther confirmed by the Chandra data presented in 
§ 2.10. 

2.6. Density, Cooling Time and Mass Pro- 
files 

2.6.1. NGC 1399 

We have derived the hot gas density radial pro- 
file of NGC 1399 from the X-ray surface bright- 
ness profile, making use of the procedure described 
by Kriss et al. (1983), that allows to deproject 
the observed emission in concentric shells, assum- 
ing spherical symmetry and homogeneity of the 
hot gas. Since the global gas distribution has a 
complex and asymmetric structure, we cannot ap- 
ply this method to the total emission; instead we 



applied it to the 'central' and 'galactic' compo- 
nents separately, whose ellipticities, derived from 
the bidimcnsional models (see previous section), 
are small. The spherical symmetry approximation 
is not valid for the cluster halo so that, in deriving 
the density profile, we assumed rotational symme- 
try around minor axis. We found that assum- 
ing symmetry around the major axis the difference 
is of the order of a few percents (except at very 
large radii where exposure correction uncertainties 
are dominant) and resulted in minor changes to 
our conclusions. The central, galactic and cluster 
components were assumed to be isothermal with 
temperatures of respectively kT=0.86, 1.1 and 1.1 
keV, in agreement with the PSPC temperature 
profiles derived by RFFJ. Buote (1999) showed 
that ASCA data suggest the presence of multi- 
phase gas within the central 5' of NGC 1399, bet- 
ter fitted by cooling flows models. However, as he 
notes, the narrower ROSAT energy range and the 
lower spectral resolution is unable to clearly dis- 
criminate between single and multi-temperature 
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Fig. 13. — Radial profiles compared to the bidi- 
mensional model profiles. Open circles represent 
HRI counts while triangles represent PSPC counts 
rescaled to match the HRI. The continuous line 
represents the best-fit model while dot-dashed, 
short-dashed and long-dashed lines represent re- 
spectively the central, galactic and cluster halo 
components. The residuals are shown in the lower 
panel in ^ units. The vertical dotted line marks 
the position of the PSPC support ring; at larger 
radii the exposure correction is poorly known. 

models (even though it reveals the presence of ra- 
dial temperature gradients) so that our assump- 
tion of isothermal components is adequate for the 
kind of analysis performed here. 

The electron density profiles of the isothermal 
components are shown in Figure 15. The den- 
sity profiles are centered on the centers of their 
respective components and thus the total density, 
shown as a thick line, does not take into account 
the centers offset. Including the displacements of 
each component with respect to the galaxy center 
results in a slightly smoother profile in the outer 
regions, due to the azimuthal averaging, shown as 
a thin continuous line. ROSAT and ASCA spec- 
tral fits (RFFJ; Buote 1999; Matsushita, Ohashi 
& Makishima 2000) revealed that additional ab- 
sorption, over the galactic value, may be required 
in the galaxy center. We thus estimated an up- 
per limit for the density of the central component 
assuming an absorption of 2.0 x 10^^ cm~^ (see 




Right Ascension (J2000) 



Fig. 14. — Residuals of the bidimensional fit on 
the 15 arcsec/pixel HRI image, after convolving 
with a gaussian oi a ~ 30". Colors from black to 
white represent residuals from — 2.5(J to 22cr. For 
comparison we superimposed the contour levels of 
the adaptively smoothed image (Figure 3a). 

RFFJ), shown as downward triangles in Figure 15. 

We must notice that adding the single density 
components, derived from the surface brightness 
decomposition, to obtain the total density profile 
is not strictly the correct approach since the sur- 
face brightness does not depend linearly on the 
gas density. However regions of the image in which 
more than one component are significant are small 
and the good agreement of our results with the 
estimate of RFFJ (who deprojected the total sur- 
face brightness) indicates that we are not introduc- 
ing significant errors in our total density estimate. 
The same considerations hold for the mass profiles 
derived below. 

Using the density profiles and the assumed tem- 
peratures we estimated the cooling time oc 
nkT /neTiH^ as a function of radius, where Ue is 
the electron density shown in Figure 15, nn = 
7ie/1.21 (Sarazin 1988) and n = rte + ^^^r arc the 
hydrogen and total density and A is the cooling 
function given by Sarazin & White (1987). Figure 
16 shows that the central and galactic components 
have cooling times much shorter than the age of 
the universe, assumed to be ~ 10^" years, out to 
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(a) (b) 

Fig. 17. — Integrated mass profiles derived from deprojection. (a) Integrated gass mass, with, symbols having 
the same meaning as in Figure 15. (b) Integrated gravitating mass for our different temperature profiles 
(see discussion in text) arc shown as continuous lines. Open squares and filled circles represent the total 
gravitating mass estimated respectively by Rangarajan et al. (1995) and Ikebe et al. (1996). The gravitating 
mass range derived by Saglia et al. (2000) from optical observations of stars, GCs and PNe, is shown as 
a shaded region. The contribution of the stellar matter measured (extrapolated) by Saglia et al. (2000) is 
represented by the thick (thin) dashed line. 



250" and 350" (23 and 32 kpc) respectively. 

The cumulative gas mass profile is shown in Fig- 
ure 17a as a function of radius. The flat profllc of 
the galactic component over 600" (55 kpc) reflects 
the steep slope of the Beta profile. The filled tri- 
angles show that the adoption of a high absorbing 
column in the galaxy center docs not affect signif- 
icantly the gas mass estimates, giving a difference 
of ~33%. We estimated the total gravitating mass 
using the expression (Fabricant, Lecar & Goren- 
stein 1980): 

, kTir) fdlnp dlnT\ 

We tried four different temperature profiles based 
on published data: we used both the linear and 
power-law approximations derived by Jones et al. 
(1997) based on PSPC spectral fits, a constant 
temperature of 1.1 keV that approximates the 
isothermal profile found by RFFJ and a power- 
law profile T{r) = Tor", with Tq = 0.6 keV and 



a = 0.13, representing the temperature drop in 
the inner 150" (14 kpc). The gravitating mass, 
extracted in concentric annuli centered on the X- 
ray centroid (i.e. taking into account the offset 
between the different components), is shown in 
Figure 17b. The similarity of the continuous lines 
indicates that the gravitating mass is weekly de- 
pendent on the assumed temperature profile and 
is mostly determined by the density profile. The 
largest deviation is seen in the model taking into 
accoimt the central temperature drop, which is 
systematically lower in the inner 100". Our pro- 
files are in good agreement with RFFJ and fall 
within the range predicted by optical observations 
(shaded region, Saglia et al. 2000). We clearly 
distinguish the transition between the galaxy and 
cluster halo near 700" (64 kpc), in agreement 
with the ASCA data (filled circles, Ikebe et al. 
1996). Both components exceed the luminous 
mass, based on an extrapolation of stellar mea- 
surements (thin dashed line), indicating the pres- 
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Fig. 15. — Density profiles derived from deprojec- 
tion. Dot-dashed, short dashed and long dashed 
lines represent respectively the central, galactic 
and cluster halo component. Downward trian- 
gles show the upper limit on the central compo- 
nent profile adopting a central absorbing column 
of 2 X 10^^ cm~^. The thick (thin) continuous 
lines show the total density ignoring (including) 
the center offset of the different components (see 
discussion in text). For comparison the values ob- 
tained from Rangarajan et al. (1995) are shown as 
filled squares. As in Figure 13 the vertical dotted 
line marks the position of the PSPC support ring. 

ence of a consistent amount of dark matter. 

Thanks to the ROSAT HRI resolution in the 
inner 100" (9 kpc) we are able to resolve the dy- 
namical behavior of the central component. The 
gravitating mass follows closely the stellar mass 
profile, confirming that the central galaxy dynam- 
ics are dominated by the luminous matter (Saglia 
et al. 2000). 

2.6.2. NGC 1404 

Following the method outlined in the previous 
section we derived the electron density, cooling 
time and mass profiles for NGC 1404. Since the 
deprojection algorithm tends to enhance the fluc- 
tuations present in the brightness profile, we used 
the best fit model to derive the deprojected elec- 
tron density. This is possible because the NGC 
1404 radial brightness profile was well fitted by a 
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Fig. 16. — Cooling time profiles derived from de- 
projection. The symbols have the same meaning 
as in Figure 15. Note that the assumption of a 
high central absorption translates into lower lim- 
its on the central cooling time (filled triangles). 

single symmetric Beta model out to 90" (§ 2.4.2). 
Figure 18a shows the good agreement between the 
deprojected density derived from the Beta model 
(continuous line) and the one from the surface 
brightness profile (filled circles). The discrepancy 
in the inner 5" is due to the fact that the Beta 
model is deconvolved for the instrumental PRF. 
We stopped the deprojection at 80" to avoid un- 
certainties related to the background subtraction 
and to the asymmetries due to ram pressure strip- 
ping (§ 2.4.2), both of which tend to increase the 
fluctuations at large radii. 

The cooling time profile derived from the Beta 
model is shown in Figure 18b and suggests that 
significative cooling must be present in the NGC 
1404 halo because the cooling time is smaller than 
10^ yr within the inner 80" (~7 kpc). We adopted 
a cooling function A = 3.8 x 10~^^ ergs cm^ s~^ 
for kT=0.6 keV (Sarazin & White 1987). 

The integrated gas mass and the total gravi- 
tating mass derived from equation 6 are shown in 
Figure 19. We assumed NGC 1404 to have an 
isothermal profile with kT=0.6 keV (Jones et al. 
1997). 
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Fig. 18. — (a) Deprojected electron density profile of NGC 1404. The filled circles represent the profile 
obtained from direct surface brightness deprojection while the continuous line shows the smoother profile 
obtained by deprojecting the best-fit Beta model, (b) Cooling time profile derived from the best fit Beta 
model. 



2.7. Global Fluxes and Luminosities 

To determine the total HRI X-ray flux of NGC 

1399 wc measured the total net counts within a 
500" radius from the X-ray centroid. The back- 
ground counts were extracted from the rescaled 
particle map (sec § 2.4.1) in the 500" 850" anmi- 
lus . We then corrected the background subtracted 
counts for the contamination of pointlike sources 
within the 500" circle. This was done by subtract- 
ing, for each source, the net counts measured by 
the wavelets algorithm (see § 2.9). The final count 
rate, for the hot halo of NGC 1399, is 0.451±0.004 
counts sec~^. 

To convert this count rate into a flux we used 
the best fit Raymond- Smith (RS) model found 
by RFFJ from the ROSAT PSPC data (a main 
thermal component with kT ~1 keV plus a softer 
emission with kT = 80.8 cV contributing for 16% 
of the emission in the 0.2-0.3 keV range). Using 
the HEASARC-XSPEC software, we simulated a 
RS spectrum fixing the absorbing column to the 



galactic value (1.3 x 10 " cm 



^) and the metal 
abundance to the solar value. With this spec- 
trum, we calculated the conversion factor between 



counts and fiuxes using the PIMMS software, that 
takes into account the HRI spectral response func- 
tion. We obtain 1 count sec~^=3.158x 10~^^ erg 
s-i cm-2 so that fx = (1.42 ± 0.01) x IQ-^^ erg 
s~^ cm^^ in the 0.1-2.4 kcV energy range. As- 
suming a distance of 19 Mpc, Ljc(0.1 — 2AkeV) = 
(5.50 ± 0.04) X 10*1 erg 8"^ Taking into account 
the uncertainties in the exposure correction and 
background extraction, this result agrees within 
10% with the estimate of RFFJ, converted to our 
adopted distance. 

In the case of NGC 1404 we used a single tem- 
perature Raymond-Smith spectrum with kT = 0.6 
keV and metal abundance 0.2 solar, following 
Jones et al. (1997). We measured a count rate of 
0.117 ±0.001 counts sec~i within a 150" radius, 
extracting the backgroimd from the 200" -300" an- 
nulus. The corresponding (0.1-2.4 keV) flux and 
luminosity are fx 
cm~^ and Lx 

We also searched for evidence of nuclear vari- 
ability of the two galaxies in our data, finding 
none. 



(4.06 ± 0.03) X 10-12 erg s" 
(1.57 ±0.01) X lO''! erg s'^. 



16 




10 1 00 

r (arcsecy 



Fig. 19. Integrated mass profiles of NGC 1404 
derived from deprojection. The total gas mass is 
represented by the lower line while the total mass 
derived using equation 6 is shown by the upper 
line. 

2.8. X-ray/Radio Comparison 

NGC 1399 is known to be a radio galaxy with a 
faint radio core (Score ^ 10 mJy at 5 GHz, Ekers 
et al. 1989). Killeen, BickneU & Ekers (1988) stud- 
ied in detail the 6 and 20 cm radio emission. The 
6 cm radio contours (their Figure 1) arc superim- 
posed on the 'csmoothed' X-ray image in Figure 
20. The radio core lies in the center of the galaxy, 
but the nuclear source is not distinguishable from 
the radio jets at this resolution. 

We searched our X-ray data for a central point 
source that could be the counterpart of the ra- 
dio core. Trying to model the nuclear emission 
with a Beta model plus a delta function, con- 
volved with the HRI PRF, resulted in minor vari- 
ations in the Beta model parameters and no im- 
provement in the fit statistics, thus indicating no 
need for additional nuclear emission. To estimate 
an upper limit we varied the central point source 
brightness and repeated the fit until we reached 
a 3(T confidence level. We obtain an upper count 
rate limit of 2.25 x 10~^ counts s~^. Assuming 
a power law spectrum with Uph = 1.7 and galac- 
tic absorption, we get fx^ = 1.0 x 10~^^ erg s~^ 



cm~^ in the 0.1-2.4 kcV band, corresponding to 
= 3.9 X 10^^ erg s~^. This value is approx- 
imately half of the one derived by Sulkanen & 
Brcgman (2001) making use of one of our HRI 
observations (RH600831a01) and 30% of the limit 
estimated by RFFJ by means of spectral analy- 
sis. Recently Locwcnstein et al. (2001) were able 
to further reduce the upper limit on the nuclear 
source by one order of magnitude making use of 
Chandra data. 

By comparing the 6 cm radio contours with the 
X-ray halo (Figure 20) we identify several inter- 
esting features; i) in the X-rays the nuclear region 
is elongated in the N-S direction, following the ra- 
dio jets. This elongation is more pronounced in 
the Southern direction, where there is an excess of 
X-ray emission following the Western side of the 
jet (region A); ii) At the point where the South- 
ern jet ends in the lobe there is a steep gradient in 
the X-ray emission. The lobe is bent towards the 
West and aligned with a region of low X-ray emis- 
sion (region B); iii) The Southern lobe ends in an 
irregular clump coincident with a region of higher 
X-ray emission (region C); iv) No features as the 
ones seen in the Southern lobe are evident in the 
Northern one. Here the radio jet ends smoothly in 
a regular lobe with no sign of sharp transition and 
there is no evidence of cavities (minima) in the 
X-ray emission. Instead the jet and the lobe are 
aligned with an X-ray structure that seems to be 
coincident with the right side of the jet/lobe. At 
the North- West end of the Northern lobe there are 
two X-ray clumps (region D): they are detected as 
a single source by the wavelet algorithm (No. 32 
in Table 3). To a more accurate visual inspection 
we found that the smaller and compact source is 
coincident with an optical (probably background) 
counterpart, while there is no counterpart for the 
diffuse clump (§ 2.9) suggesting that it is a local 
feature of the hot gaseous halo. 

From the residual image of the two-dimensional 
fit (Figure 14) we can calculate the significance of 
these features. We find that the depiction corre- 
sponding to the Southern lobe (region B) and the 
Northern diffuse clump (region D) are significant 
respectively at the 2.5a and ~ 4(7 level. Region 
C corresponds to a 3a brightness enhancement (§ 
2.5). 

To examine in greater detail the nuclear region 
we subtracted the halo model developed in § 2.5 



17 



3:38:30 



20 



Right Ascension (J2000) 



Fig. 20. — 6 cm radio contours (green; from Figure 1 of Killeen, Bicknell & Ekers 1988) superimposed on the 
adaptivcly smoothed 1" /pixel X-ray image. The logarithmic grayscale represents intensities from 8.6xl0~^ 
(black) to 6.7x10"^ (white) cnts arcmin"^ s~^. 



from the X-ray image to produce a residual map 
of the nuclear region with a 2" /pixel resolution 
(Figure 21). As discussed at the beginning of this 
section, there is no evidence of a nuclear point 
source aligned with the radio centroid. Instead 
residual emission is visible on both sides of the 
radio jets a few arcseconds (< 1 kpc) from the 
nucleus. There is no evidence that any of these 
excesses may be due to the nuclear source becaiisc 
they are displaced with respect to the radio emis- 
sion. Moreover the radio emission seems to avoid 
the X-ray excess on the eastern side of the nucleus, 
suggesting that these features may be correlated. 

In coincidence with the direction of the radio 

jets there are two regions of negative residuals 
(~ —2a), extending up to ~ 20", where additional 
excess emission, on both the Northern and South- 
ern sides of the X-ray core and aligned with the 
radio jets, is clearly visible. The Southern excess 



corresponds to region A of Figure 20. Residual 
emission is also present on the Western side of the 
nuclear region. Visual inspection of the Chandra 
image (see § 2.10) shows that this excess is proba- 
bly due to the presence of point sources unresolved 
in the HRI image. 

From the density profiles obtained in § 2.6.1 
we were able to obtain the pressure profiles of 
the three halo components shown in Figure 22. 
The pressure profiles are compared to the mini- 
mum radio pressure measured by Killeen, Bick- 
nell & Ekers (1988) in the Northern (circles) and 
Southern (triangles) radio jets/lobes. The Figure 
confirms their result that the radio source can be 
confined by the ISM whose pressure is higher than 
the radio one. We must notice that they found a 
difference between radio and thermal pressure of 
more than one order of magnitude, while we find 
Pt hernial /PRadio = 3-=-4. This is mostly due to the 
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3:38:30 27 

Right Ascension (J2000) 

Fig. 21. — 6 cm radio contours superimposed on 

the 2" /pixel residuals X-ray image of the nuclear 
region. The image has been smoothed with a gaus- 
sian of fT = 4". The logarithmic grayscale repre- 
sents residuals ranging from —2.2a (black) to +4a 
(white). The circles on the left show the positions 
of X-ray sources detected by the wavelets algo- 
rithm (No. 2 and 11 in Table 3). The radius of the 
central source (corresponding to the central X-ray 
peak) is ~ 22". 

poor ISM temperature constraint obtained from 
the Einstein IPC data and to our better spatial 
resolution. 

The vertical error bar in the Figure indicates 

the pressure estimate for the diffuse clump ob- 
served in the D region. To calculate this value 
we assumed that the brightness enhancement is 
due to a homogeneous sphere of ^ 9" (800 pc) ra- 
dius at the same temperature of the galactic halo 
(1.1 keV). Depending on whether the background 
is extracted locally (an annulus 9" wide) or from 
the Beta model, the 'bubble' pressure can vary 
up to 25%; so that the allowed range is shown as 
an error bar. This estimate is dependent on the 
assumed geometry of the emitting region and on 
projection effects, so that it must be considered as 
an upper limit on the actual pressure. 



10 
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Fig. 22. — Gas pressure profiles derived from de- 
projection. The symbols have the same meaning of 
Figure 15. Open circles and triangles represent re- 
spectively the radio pressure of the Northern and 
Southern lobes measured by Killeen, Bicknell & 
Ekers (1988). The vertical error bar in the Figure 
indicates the pressure range of the D region (see 
discussion in text). 

2.9. Discrete Sources 

We used the algorithm developed at the Palermo 
Observatory by F. Damiani and collaborators 
(Damiani et al. 1997a) to detect the discrete 
sources present in the HRI field. The algorithm 
measures the relevant parameters in a wavelet 
transformed space, using different wavelet scales 
to optimize the detection of both compact and ex- 
tended sources. The wavelet algorithm found 43 
sources in our field, using a signal to noise thresh- 
old of 4.65. The threshold was chosen to assure 
a contamination of one/two spurious detections 
per field (F. Damiani private communication, see 
also Damiani et al. 1997b). Three more sources 
were found with the algorithm, but being clearly 
associated to the ribbon-like features due to the 
presence of hot spots in the detector (Figure 2), 
were excluded from the subsequent analysis. Af- 
ter visual inspection of the adaptively smoothed 
image (Figure 3) three more detections (No.44, 
45, 46) with a lower S/N ratio but showing clear 
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Right Ascension (J2000) 

Fig. 23. — The sources detected with the wavelets algorithm of Damiani et al. (1997a, b) superimposed on 

the 5'' /pixel adaptivc^ly smoothed HRI image. Solid circles show the region of maximum S/N ratio for each 
source. Dashed circles represent reff and 5re// for NGC 1399 and NGC 1404 (Goudfrooij et al. 1994). 



features of pointlike sources, were added to the 
list. 

The source positions in the field are shown in 
Figure 23. Sources No. 42 and 43 are clearly ex- 
tended and seem to be associated to structures 
in the NGC 1399 halo revealed by the adaptive 
smoothing algorithm (§ 2.3). 

Table 3 lists the parameters measured with the 
wavelets algorithm for each source: position, 3ct 
radius, number of counts with statistical error, 
maximum signal to noise ratio (for sources de- 
tected at different spatial scales), count rate and 
flux in the 0.5-2.0 keV band for both HRI and 
PSPC data. The 0.5-2.0 keV energy band was 
chosen to allow an easily comparison with the 
Hasinger et al. (1993) background counts (see be- 
low). The the last three columns of Table 3 report 
the source short-term variability obtained with a 



Kolmogorov-Smirnov test of the HRI light curves, 
long-term variability from comparison with PSPC 
fluxes (as explained below), and possible opti- 
cal counterparts (within the 3a radius measured 
by the wavelets algorithm) found by both visual 
inspection of the DSS photographic plates and 
cross-checking with known sources present in the 
Nasa Extragalactic Database (NED). The names 
of the counterparts correspond to those reported 
by Hilker et al. (1999a) or Smith et al. (1996) in 
their optical and NUV studies of the Fornax field. 
Sources No. 27 and 28, representing respectively 
NGC 1399 and NGC 1404, are not included in 
the table. 
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Table 3: Discrete Sources Properties 



Source 
No. 


R.A. Dec. 
(J2000) 


3ct Radius 
(arcsec) 


Counts 


Max S/N 


Count Rate 
(_1U Onts sec J 


m (0.5 2.0 keV) 

/in — 15 ~ — 1 — 2\ 
(^iU erg s cm ) 


fmax abs. 

H RI (0.5 — 2.0 keV) 

/in — 15 _ — 1 „. — 2\ 

(iU erg s cm j 


f.pOlU . i ft If 

P S PC{0.5 2.0 keV) 

fin — 15 _ — 1 — 2\ 
(lU erg s cm ) 


Time Var.* 
snort-term / long-term 


Opt. id. 


NED obj.^ 


1 


3:38:20.0 


-35:24:42 


8.7 


35±12 


4.7 


2. 1±0.7 


4.3±0.2 


5.6 


28±7 


n/y 






2 


3:38:32.0 


-35:26:45 


8.7 


54±17 


6.1 


3.2±1.0 


6.7±0.3 


8.6 




n/- ■ ■ 






3 


3:38:25.3 


-35:22:42 


8.7 


49it 14 


6.7 


2.9±0.8 


6.1±0.3 


7.8 




n/- ■ ■ 






4 


3:38:45.9 


-35:22:50 


8.7 


42ilz 12 


6. 1 


2.6±0.7 


5.3±0.3 


6.9 


9±3 


n/n 


y 


CGF 0202 




3:38:31.8 


-35:26:02 


11.9 


165 ib23 


13.0 


10±2 


20±1 


26.0 




n/- ■ ■ 




CGF 0233 




3:39:12.3 


-35:28:12 


13.2 


172±21 


14.9 


10±1 


21±1 


27.5 


25±4 


n/n 


y 




7 


3:38:31.8 


-35:33:42 


12.2 


o6±16 


6.0 


3.3±0.9 


6.9±0.3 


8.9 


8±3 








8 


3:38:41.2 


-35:31:35 


10.8 


1000ib35 


60.6 


59±2 


123±6 


159.2 


170±8 


99%/y 


y 




9 


3:38:49,8 


-35:24:57 


12.2 


47±14 


5.0 


2.8±0.9 


5.8±0.3 


7.5 




n/- ■ ■ 


" 




10 


3:39:08,1 


-35:31:22 


12.2 


o6±lb 


6.2 


3.3±0.9 


6.9±0.4 


8.9 


22±6 


n/n 






11 


3:38:32.4 


-35:27:05 


13.7 




11.5 


10±2 


21 ±1 


27.2 




n/- ■ ■ 






12 


3:38:51.7 


-35:26:42 


11.6 


295 zh22 


24.7 


r _^ 


36±2 


47.0 


58±6 


95%/y 


y 


yjyjrP UlUz 


13 


3:38:48.6 


-35:28:35 


12,2 


'^''^^^^ 


1 


4.5±1.1 


9.3±0.5 


12.0 




n/- ■ ■ 


y 




14 


3:39:09.7 


-35:27:07 


12.2 


46 ±14 


5 . 


2.8±0.8 


5. 7±0.3 


7.4 




n/- ■ ■ 


y 




15 


3:38:01,9 


-35:26:27 


17,3 




4.7 


3.6±1.1 


7.5±0.4 


9.7 




n/- - ■ 






16 


3:38:29.4 


-35:25:05 


17.3 


75 ±23 


4.8 


4.4± 1.3 


9. 1±0.5 


11.7 




n/- ■ ■ 


" 




17 


3:38:31.6 


-35:30:58 


12.2 


70±18 


6.8 


4.1±1.1 


8.6±0.5 


11.1 




n/- ■ ■ 


" 




18 


3:38:03.7 


-35:33:27 


17.3 


65± 19 


5.4 


4.1±1.2 


8.5±0.4 


10.9 




n/- ■ ■ 


y 


CGF 0354 


19 


3:38:08.5 


-35:34:20 


15.0 


908 ±34 


52.3 


57±2 


119±6 


152.7 


127±7 


99%/n 


y 




20 


3:38:43.5 


-35:33:50 


15.2 


105± 18 


8.7 


6.3± 1. 1 


13. 1±0. 7 


16.8 




n/- - - 


y 






3:39:20.1 


-35:29:22 








32±2 


68±4 


86.9 


83±7 


99%/n 


y 




22 


3:38:24.1 


-35:14:52 


17.3 


60± 18 


4.9 


3.8± 1.2 


8.0±0.4 


10.3 


12±4 


n/n 






23 


3:38:12.2 


-35:16:17 


17.3 


59±18 


4.8 


3.7±1.1 


7.7±0.4 


9.9 




n/- - - 


y 




24 


03:37:58.5 


-35:19:42 


16.7 


379±27 


24.9 


24±2 


48±2 


61.6 


75±7 


n/y 


y 




25 


03:37:34.8 


-35:24:04 


17.3 


64±19 


5.3 


4.0±1.2 


8.3±0.4 


10.6 




n/- - . 






26 


03:38:43.1 


-35:33:10 


17.3 


73±21 


5.6 


4,3±1.2 


9.0±0.5 


11.6 




n/- . - 


y 




29 


03:39:43.0 


-35:24:11 


24.3 


143±27 


8.6 


10±2 


21±1 


26.7 


62±8 


n/y 


y 




30 


03:38:31.0 


-35:14:17 


26.4 


586±38 


26.5 


38±2 


79±4 


102.1 




99%/- ■ - 


y 


FCCB 1263 


31 


03:38:17.1 


-35:15:57 


24.5 


124±31 


7.2 


8±2 


16.4±0.8 


21.1 




n/- ■ ■ 


y 




32 


03:38:25.7 


-35:25:27 


24.5 


145±37 


6.2 


8.5±2.2 


17.7±0.9 


22.8 


27±7 


n/n 


y 




33 


03:37:36.5 


-35:33:37 


24.2 


921±40 


41.4 


60±3 


125±6 


161.5 




n/... 


y 




34 


03:37:39.4 


-35:34:17 


21.9 


426±31 


24.2 


28±2 


58±3 


74.5 




n/... 


y 




35 


03:37:50.1 


-35:21:27 


34.6 


139±37 


5.8 


8.4±2.3 


17.6±0.9 


22.6 




n/... 


y 




36 


03:39:27.8 


-35:18:52 


34.6 


129±35 


5.6 


8.6±2.3 


18.0±0.9 


23.2 


46±8 


n/y 


y 




37 


03:37:36.3 


-35:14:47 


49.0 


153±44 


5.1 


11±3 


24±1 


30.7 


25±7 


n/n 






38 


03:37:33.4 


-35:17:17 


49.0 


190±49 


6.0 


13±3 


27±2 


35.2 




n/... 






39 


03:37:29.0 


-35:41:11 


69.3 


580±119 


14.1 


67±14 


140±7 


180.1 


232±12 


n/... 


y 




40 


03:38:14.7 


-35:17:57 


69.3 


244±67 


5.2 


15±4 


31±2 


40.3 




n/... 






41 


03:38:36.7 


-35:41:32 


62.9 


366±61 


8.9 


24±4 


51 ±3 


65.8 


57±7 


n/n 






42 


03:39:12.3 


-35:26:12 


98.0 


312±91 


4.7 


19±5 


39±2 


50.7 




n/... 






43 


03:38:12.2 


-35:30:02 


98.0 


365±101 


5.0 


22±6 


46±2 


59.1 


60±14 


n/n 






44 


03:38:05.1 


-35:19:05 


12.2 


42±13 


4.6 


2.5±0.8 


5.3±0-3 


6.8 




n/... 






45 


03:37:59.9 


-35:23:30 


17.3 


56±18 


4.4 


3 . 3 ± 1 . 


6.8±0.4 


8.8 




n/... 






46 


03:37:59.2 


-35:11:37 


49.0 


134±42 


4.4 


9, 9 ±3,1 


20±1 


26.6 


28±8 


n/n 







"Short-term = probability of the Kolmogorov-Smirnov test: 'n' means no variability detected at the 95% level. 
Long-term = difference between PSPC and HRI counts: 'y' means variable source (> 3cr), 'n' means no variability detected, '• • •' means no available 

data. 

''CGF=Hilker et al. (1999a); FCCB=Smith et al. (1996) 



We tried to determine which X-ray sources be- 
long to the Fornax cluster both by visual inspec- 
tion of the photographic DSS plates and by com- 
paring their positions with those of published cat- 
alogs of the Fornax region. Sources No. 2 and 
11 are likely to belong to NGC 1399 because 
of their position within the optical effective ra- 
dius (re// = 44.7", Goudfrooij et al. 1994) while 
sources No. 16 and 32 appear to be coincident 
with the position of two globular clusters (Kissler- 
Patig et al. 1999). Two sources are clearly as- 
sociated with background objects: source No. 4 
corresponds to a background galaxy located at 
z = 0.1126 and represents the X-ray counterpart 
of the strongest component of the radio source 
PKS 0336-35 (Carter & Malin 1983); source No.30 
is coincident with two optical background galax- 
ies classified as interacting by Ferguson (1989) and 
clearly visible on the DSS plates. For the remain- 
ing sources there is no clear indication of whether 
they belong to the Fornax cluster or not. 

Sources No. 8, 12, 19, 21, 30 are variable at the 
95% confidence level, according to our KS tests. 
Source No.30 is associated with a short X-ray 
burst occurred in the background galaxy between 
January and August 1996. In fact the source 
is not detected in both the August 1991 PSPC 
(RP600043n00) and the February 1993 HRI ob- 
servations (RH600256n00). It is present, instead, 
in the last two HRI observations (RH600831n00 
and RH600831a01) with a decrease in the observed 
count rate of a factor 8 between the two. 

We must notice that our background level was 
also variable so that the presence of temporal vari- 
ability was clearly detect only for the strongest 
sources (/x(o.5-2.o keV) > 'Sxli)-^^ ergs^i cm^^)^ 
However, since the background subtraction was 
performed locally (i.e. extracting the background 
in a region surrounding the source, see Damiani 
et al. 1997a), this variability does not affect the 
source counts reported in Table 3. To further 
check the reliability of the wavelets algorithm, we 
extracted the counts manually for each source us- 
ing a local annulus to measure the background, 
finding consistent results (within la) with those 
derived by the wavelets technique. 

We cross-checked our source list with the list 
of sources detected on the PSPC image by the 
GALPIPE project (see §2.5). Several sources de- 
tected in the HRI field were bright enough and 



sufficiently isolated to allow a comparison of the 
PSPC flux with the HRI results and to check for 
long-term variability. We converted the PSPC 
count rates to fluxes in the 0.5-2.0 keV band (1 
cps=6. 397x10"^^ ergs cm~^ s~^), assuming the 
same power law spectral model used for convert- 
ing the HRI counts (sec below). The PSPC fluxes 
and long-term variability (i.e. when PSPC and 
HRI flux estimates differ more than 3f7) are in- 
cluded in Table 3. 

This comparison also confirmed that sources 
No. 39 and 46, both near the very edge of the HRI 
FOV, are not an artifact of the poor S/N ratio 
or the steep emission gradient because both are 
present in the PSPC data as well and detected 
with high significance. 

We compared our source flux distribution to 
the expected background counts obtained from a 
deep survey of the Lockmann hole (Hasinger et al. 
1993). We used the sensitivity map generated by 
the wavelet algorithm to weight each source for 
the effective detector area above the correspond- 
ing limiting flux. In this way we are able to nor- 
malize the measured counts at each flux to the to- 
tal surveyed area. We computed the flux of each 
source assuming a power law spectrum with pho- 
ton index 1.96±0.11 {Jh°riIo%-2.o keV) Table 3), 
as found by Hasinger and collaborators for faint 
sources. This gives a conversion factor of 1 cut 
s-^=2.087 X 10"" ergs s"^ cm'^ in the 0.5-2.0 
keV band (the same used by Hasinger). We also 
tried a Raymond-Smith spectrum with kT = 0.52 
and abundance ~ 0.2 solar, following the results of 
Kim & Fabbiano (1995) in their study of pointlike 
sources near NGC 507 and obtaining an almost 
identical result. 

The flux distribution is shown in flgure 24 where 
our differential counts (empty circles) are com- 
pared to the Hasinger ones (continuous line). The 
horizontal error bars take into account uncertain- 
ties in the the power law slope. Arrows represent 
the upper limit obtained flxing the absorption to 
the upper limit of 4.5 x 10^° cm'^ (/^"(^"g^^a'o keV) 
in Table 3), allowed from the spectral analysis 
performed by RFFJ. Our counts are higher than 
the Hasinger estimate, giving a total of 312 ± 47 
sources vs 230 ± 15 expected counts, but still com- 
patible within 2a. To take into account absorption 
uncertainties we estimated a lower limit on the 
number of expected sources adopting the higher 
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(a) (b) 

Fig. 24. — (a) Measured vs expected source counts. Measured counts, corrected for sensitivity variations 
across the detector, are represented by empty circles, while the continuous line is the expected distribution 
from Hasinger et al. (1993). Horizontal error bars take into account uncertainties in the power law slope. 
Arrows represent upper limits obtained fixing the absorption to 4.5 x 10^° cm~^ (Rangarajan et al. 1995). (b) 
Radial excess of sources over the expected number from Hasinger et al. (1993). The excess is concentrated 
around NGC 1399, thus suggesting to be due to sources associated with the galaxy. 



absorption value found by RFFJ in the NGC 1399 
field. In this case the number of expected counts 
decreases to 184 ± 14 leading to an upper limit of 
the observed excess of ~ 2.6a. 

The spatial distribution of detected sources is 
not uniform. In Figure 24 we show the source 
excess in circular annuli centered on NGC 1399. 
We find that the excess is peaked on the domi- 
nant galaxy, suggesting that the inner sources are 
likely to be associated with NGC 1399 rather than 
being background objects. In fact, excluding the 
central 300" (28 kpc), the measured counts drop 
to 244 ± 37, in agreement within less than la 
from the expected value. The X-ray luminosity 
of the sources within the central 300" range from 
1.7 X 10^*^ erg s^^ to 1.4 x 10^^ erg s"\ if they are 
in NGC 1399. Thus all of them have luminosities 
in excess of 1.3 x lO"*^ erg s^^, the Eddinghton 
luminosity for a solar mass accreting object, sug- 
gesting that they may have massive Black Hole 
companions if they are accretion binaries. 



2.10. The Chandra data 

While this work was in progress the Chandra 
data on NGC 1399 became public. A fidl data 
reduction of the Chandra data, which is how- 
ever already in progress (Angelini, Loewenstein & 
Mushotzky 2001; Loewenstein et al. 2001), is be- 
yond the aim of the present work. However we de- 
cided to compare the ROSAT data with the Chan- 
dra image to see if our conclusions are supported 
by the improved Chandra resolution and sensitiv- 
ity 

We used the longest Chandra observation of 60 
ks (Obs. ID 319, 55.942 ks live time), performed 
on 2000 January 18 with the ACIS-S detector. We 
produced a 0.5 arcsec/pixel image in the energy 
band 0.3-10 keV, from the ACIS-S3 chip, which 
covers the central 8x8 arcmin of NGC 1399. 
The image was then adaptively smoothed with the 
csmooth algorithm included in the CIAO package. 
The final image is shown in Figure 25. No expo- 
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Fig. 25. — Adaptively smoothed Chandra ACIS- 
S image of NGC 1399 in the energy range 0.3-10 
keV. Colors from black to white represent count 
rates from 0.021 cnts s^^ arcmin"^ to 7.7 cnts s^^ 
arcmin^^. Circles represent the position of sources 
detected by the wavelets algorithm in the HRI im- 
age. 

sure correction was applied because the csmooth 
algorithm works better on the uncorrected data. 
However the exposure map is quite flat and we 
further checked that any feature present in the ex- 
posure map does not affect our conclusions. 

Figure 25 shows a large number of pointlike 
sources present in the NGC 1399 field, most of 
which were unresolved in the HRI image. We 
marked with circles the sources detected in the 
HRI field (sec § 2.9). The Chandra data show that 
sources 2, 11 and 5, identified as single sources in 
the HRI image, are instead multiple sources. 

The filamentary structures found in the HRI 
image (Figure 4) are present in the Chandra image 
as well. The 'arm' protruding on the Western side 
and the 'voids' in the X-ray emission are clearly 
seen. The hypothesis that the filamentary struc- 
tures may be due to the presence of a large number 
of pointlike sources unresolved in the HRI image, 
is ruled out by Chandra data. In fact the csmooth 
algorithm clearly separates point sources from the 
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Fig. 26. — Raw 4 arcsec/pixel Chandra image of 
NGC 1399. The long arc following the Western 
edge of the radio lobe is visible South of NGC 
1399. 

extended structures (Figure 25) and shows that 
the two components are not spatially correlated. 

The correspondence between the halo features 
and the radio jets/lobes, discussed in § 2.8, are 
confirmed. Moreover the Chandra data revealed 
the presence of a long arc extending South that 
follows closely the Western edge of the radio lobe. 
The arc is better visible on the raw Chandra image 
shown in Figure 26. This feature has a low statis- 
tical significance (2(j) but the correlation with the 
radio lobe argues in favor of a real structure. 

3. DISCUSSION 

The study of galactic and cluster gaseous ha- 
los has often been based on the assumption of 
homogeneity and spherical symmetry. The deep 
HRI image of the NGC 1399 halo has shown in- 
stead significant departures from these simplified 
assumptions. Here we briefiy summarize the result 
of our analysis, and then discuss their physical im- 
plications. 

We model the X-ray brightness distribution of 
NGC 1399 with three distinct components: (1) a 
central component (dominating for r < 50"), co- 
incident with the optical galaxy; (2) a 'galactic' 
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component (50" < r < 400") displaced South- 
West of NGC 1399 and (3) an ellipsoidal 'clus- 
ter' component (r > 400") centered North-East 
of the optical galaxy. The cooling times of the 
central and galactic components are smaller than 
the Hubble time, suggesting the presence of a cen- 
tral cooling flow. The binding mass estimates, ob- 
tained assuming hydrostatic equilibrium, revealed 
that the dynamics of the central component are 
dominated by stellar matter while, at radii larger 
than 2 arcmin there is a large dark matter halo 
with two distinct components dominating respec- 
tively on galactic and cluster scales. 

We also found significant filamentary structures 
and cavities in the X-ray emission. Some of these 
structures are correlated with the morphology and 
the position of the radio jets and lobes of NGC 
1399, suggesting interactions between the radio 
emitting plasma and the hot gas. 

For what concerns NGC 1404 the surface 
brightness distribution is well represented by a 
single Beta profile out to 90" . At larger radii the 
effect of ram pressure stripping becomes evident 
as a sharp brightness cutoff in the North- West 
sector and an elongated tail on the opposite side 
of the galaxy. The cooling times derived by the 
deprojected density profile are shorter than the 
Hubble time inside the NGC 1404 halo. The mass 
profile shows no evidence of multiple components. 

Finally we used a wavelets based algorithm to 
detect 46 sources in the HRI field, 8 of which show 
significant time variability. Their distribution is 
consistent with the background source distribution 
at the 3(7 level. However there is excess over back- 
ground counts centered on NGC 1399. The pres- 
ence of a large number of pointlike sources near 
the galaxy center is further confirmed by Chandra 
data. 

3.1. Origin of the Inner Peak of the X-ray 
Surface Brightness (r < 50") 

The center of the gaseous halo is occupied by 
a strong emission peak (see Figure 3a) centered 
on the optical galaxy. Such emission is often 
found in cooling fiow galaxies (Canizares, Stewart 
& Fabian 1983; Kim & Fabbiano 1995; see also 
Fabian, Nulsen & Canizares 1991) because of the 
higher density of cooler hot gas in the galactic cen- 
ter. RFFJ revealed the presence of a cooler region 



within the central 1.5 arcmin (14 kpc); this is com- 
patible with the extent of the region dominated by 
the X-ray peak, as shown by the radial profile in 
Figure 13. Moreover both the central and galactic 
components have cooling times below the Hubble 
time (assumed to be ~ 10^° years) within respec- 
tively 250" (23 kpc) and 350" (32 kpc), so that 
significant cooling must have taken place. The 
cooling time of the cluster component, instead, is 
compatible with the Hubble time, suggesting that 
the cooling fiow is mainly of galactic origin. 

The alternative hypothesis that the central 
peak can be due to the presence of a nuclear source 
(Jones et al. 1997) is clearly ruled out by our data, 
because the observed profile (Figure 5) is wider 
than the HRI PRF even taking into account the 
PRF uncertainties. The difference between the 
optical (with slope a ~ —1.6, Hilker et al. 1999b) 
and X-ray profiles (a ~ — 2 for the central compo- 
nent, § 2.4.1) seen in Figure 9 also suggests that 
is unlikely that the central peak could be due to 
stellar sources. Allen, di Matteo & Fabian (2000) 
analyzing ASCA data, found that the addition of a 
power-law component to their NGC 1399 spectral 
model significantly improves the quality of the fit. 
They notice that the luminosity of the power-law 
component is compatible with the contribution 
expected from a population of unresolved hard 
X-ray binaries. However, the HRI count rate that 
would be produced by this power-law component 
falls more than an order of magnitude below the 
inner component count rate. We can thus exclude 
that this hard discrete population dominates the 
central component emission. 

The gravitational mass estimate derived from 
the hydrostatic equilibrium (equation 6) shows 
that the hot gas within 100" follows closely the 
stellar mass profile. This result is independent on 
the assumed temperature profile (Figure 17) being 
mostly determined by the hot gas density profile. 
These X-ray data are thus in excellent agreement 
with the optical mass determinations of Saglia et 
al. (2000) who find no evidence of dark matter in 
the galaxy core. 

The origin of the central density enhancement 
may be due to hot gas inflow within the inner 1 
arcmin. The deprojection technique described in 
§ 2.6.1 allowed to derive the cumulative gas mass 
profile of the central component (Figure 17). Ac- 
cording to the mass deposition rate measured by 
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RFFJ, which is in good agreement with our data, 
we expect to find ~ 2.5 x 10^ M© of deposited 
gas within the central 60" (~ 5 kpc) over 10^° 
yr, while we only observe ^ 2 x 10^ M0 in the 
hot phase. Taking into account a possible larger 
than line-of-sight central absorption, the observed 
mass would increase of ^ 30%, but the deposition 
rate would also increase (see RFFJ), so that the 
result remains almost unchanged. This difference 
suggests that ^ 90% of the hot gaseous mass that 
reached the inner regions of NGC 1399 has cooled 
enough to become invisible in X-rays. However 
RFFJ argue that the missing mass is not seen in 
any other phase. The situation gets even worse 
considering that we are ignoring stellar mass losses 
which must be considerable since the metallicity of 
the ISM in the galactic center is consistent with 
stellar values (Buote 1999; Matsushita, Ohashi & 
Makishima 2000). If the age of the cooling flow 
is lower than 10^° yr, this disagreement would be 
reduced. Equating the observed hot gas mass to 
the amount expected from the cooling flow, we 
estimate a lower limit of 10^ years to the age of 
the cooling flow. However, we must notice that 
the latest studies of cooling flows in clusters and 
galaxies with Newton-XMM are showing that the 
classical cooling flow model is not able to repro- 
duce the observed temperature distributions (Pe- 
terson et al. 2001; Kaastra et al. 2001; Tamura et 
al. 2001). In this case no conclusion can be drawn 
from our data until the cooling flow physics will 
be better understood. 

3.2. Dynamical Status of the Halo 

3.2.1. Gravitational Effects 

The X-ray brightness profiles (Figure 9) shows 
that past 1 arcmin the hot gas distribution is more 
extended than the stellar distribution. If the gas 
is in hydrostatic equilibrium, additional matter 
is required to explain this broader distribution. 
The change in slope of the X-ray brightness profile 
for r > 400" further requires the presence of two 
components that we identified with the 'galactic' 
and 'cluster' halo. This complex distribution was 
previously foimd by Ikebe et al. (1996) using the 
ASCA data, and explained in terms of a differ- 
ent concentration of dark matter on galactic and 
cluster scales (Ikebe et al. 1996; Makishima et al. 
2001). Our data further show that this galactic 



component is not correlated with stellar matter, 
which is even more concentrated (Figure 17b) and 
dominates the dynamics of the galactic core. 

The kinematics inferred from the GCs and PNe 
(see Figure 6 of Kissler-Patig et al. 1999) indi- 
cate that for 1' < r < 16' (5 kpc< r <90 kpc) 
from the galaxy center, the velocity dispersion in- 
creases with radius. This result is consistent with 
our measurements and suggests either the pres- 
ence of a significant amount of dark matter asso- 
ciated with the outer radii of NGC 1399 or equiv- 
alently that the outer galactic dynamics is domi- 
nated by the cluster potential. 

We must notice however that while on galactic 
scales the mass profiles may be explained by the 
presence of dark matter, the deviations from sym- 
metry, i.e. the offset between the stellar compo- 
nent (consistent with the central component) and 
the galactic dark halo, are hardly explained in the 
galaxy formation framework. For this reason in 
next sections we take into account possible depar- 
tures from equilibrium. 

3.2.2. Galaxy-Galaxy Interactions 

An alternative scenario has recently been ad- 
vanced by Napolitano, Arnaboldi & Capaccioli 
(2001, in preparation), that could explain our 
results by means of galaxy interactions. Their 
calculations show that the energy injection due 
to galaxy-galaxy encounters can account for the 
outer galaxy dynamics as seen in PNe and GCs, 
and is in agreement with the observed X-ray tem- 
perature profile. In this case there would be no 
need for additional dark matter in the 'galactic' 
halo. Support for the tidal interaction picture 
comes from Kissler-Patig et al. (1999) who showed 
that the GC system of NGC 1399 is 10 times more 
abundant than those of neighbors galaxies, while 
their properties are very similar. In the hypothesis 
that all galaxies initially had the same number of 
GCs per unit luminosity, the excess around NGC 
1399 may be compatible with tidal stripping of 
GCs from nearby Fornax galaxies. 

The effects of tidal interactions on the gaseous 
halos of elliptical galaxies have been explored 
by D'Ercole, Recchi & Ciotti (2000, hereafter 
DRC). They found that the interaction strongly 
affects the gas flow, altering the ISM density and 
thus increasing the cooling speed and luminos- 
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ity. This mechanism is enhanced when energy in- 
jection from SNe la is included and the gas re- 
mains turbulent independently from the epoch of 
the encounter. In this process the surface bright- 
ness is scarcely affected and the global isophotes 
remain circular, although locally distorted. The 
DRC model predicts that the gas inflow produces 
steeper profiles in the central regions of the galaxy 
in the moment of maximum luminosity, indicating 
that the tidal perturbations may contribute to the 
formation of central brightness peaks such as the 
one seen in the NGC 1399 core. The steep slope of 
the galactic component (Table 2) too may be due 
to tidal gas stripping because the galaxy under- 
goes significative mass losses (up to 50%) in the 
outer regions. However, we must notice that in 
the case of NGC 1399, the effect of tidal stripping 
may be smaller than what expected from the DCR 
models because PSPC data (Jones et al. 1997) sug- 
gest a lower SN rate than assumed by D'Ercole and 
collaborators. 

For these models to be consistent with the X- 
ray observations, tidal interaction must strongly 
affect the halo density distribution: the increase 
in the M/L ratio past 2 arcmin (Figure 17) is 
mainly determined by the shallower slope of the 
gas density profile in the outer galactic regions, 
since the isothermal model still exhibits the multi- 
component behavior (§ 2.6.1). Thus a simple tem- 
perature increase would not produce the observed 
agreement with the optical data. Such strong pro- 
file alterations are not seen in the DRC simula- 
tions. Conversely Barnes (2000) found that en- 
counters between galaxies of different masses can 
produce both radial profiles with a central peak 
and a 'shoulder' similar to our multi-component 
structure and azimuthal asymmetries as those seen 
in the NGC 1399 galactic halo. However his sim- 
ulations also show that the stellar profile follows 
the gas profile, which is not seen in our data. 

The X-ray surface brightness analysis (§ 2.3) 
revealed the presence of complex structures in the 
galactic halo. We have shown in § 2.5 that these 
brightness fluctuations are not due to artifacts of 
the smoothing algorithm, because they are present 
in the residuals of the bidimensional model, and 
have high statistical significance. Due to the weak 
dependence of the HRI spectral response on the 
temperature around 1 keV (see Clarke, Harris & 
Carilli 1997), they must be due to local density 



fluctuations of the hot gas. While some of these 
features can be due to the interaction between the 
radio jets and the surrounding ISM (§ 2.8), most 
of them are not related to any evident radio or op- 
tical feature and are spread all around the galaxy, 
suggesting a 'disturbed' nature of the whole halo. 

It is possible that these structures are the sig- 
nature of interaction between the NGC 1399 halo 
and other Fornax galaxies. In fact DRC predict 
significative density fluctuations in the hot gas and 
the formation of filamentary structures near the 
galaxy center as a consequence of galaxy-galaxy 
encounters. It's not easy to determine if these 
perturbations are due to the nearby NGC 1404 
or to more distant members because, according to 
DRC, the perturbing effects lasts several Gyrs. We 
searched for possible signatures of interactions be- 
tween NGC 1399 and NGC 1404 looking for excess 
emission in the region between the two galaxies. 
The slight excess we found is only significant at 
the 2a level, so that no significant conclusion can 
be drawn from our data. 

3.2.3. Ram Pressure 

A different scenario that could explain the clus- 
ter and galactic halo displacement is one in which 
NGC 1399 is moving in the Fornax cluster poten- 
tial, thus experiencing ram pressure from the clus- 
ter hot gas. The radial profiles in Figure 9 show a 
high degree of asymmetry with respect to the opti- 
cal galaxy. The 2D models indicate that the galac- 
tic component is displaced ^ 1' (5 kpc) South- 
West with respect to the X-ray peak (and thus op- 
tical galaxy), while the cluster component seems 
located ^5.6 arcmin (31 kpc) North-East of NGC 
1399. Jones et al. (1997) argue that this displace- 
ment is compatible with the hypothesis that the 
galaxy is slowly moving (40 km s"-'^) South-East 
in a larger potential. This would explain the dis- 
placement of the galaxy with respect to the cluster 
component but not the one of the galactic halo. In 
the ram pressure scenario, the displacement and 
large scale asymmetry of the 'galactic' halo may 
suggest that NGC 1399 is not sitting at the very 
center of the cluster potential and is slowly moving 
North-East toward the 'cluster' component. 

The effect of ram stripping has already been 
observed in nearby galaxies. For instance Irwin & 
Sarazin (1996) observed the presence of a very ir- 
regular halo in the Virgo galaxy NGC 4472. They 
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found a steeper gradient in the X-ray emission in 
the direction of the cluster center associated with 
a temperature enhancement. They calculate that 
ram stripping is able to explain the observed halo 
features if the galaxy is moving at ~ 1300 km 
in the cluster potential. Their ROSAT HRI images 
also revealed filamentary structures and holes in 
the galactic halo, very similar to the features that 
we observed in the adaptively smoothed image of 
the inner halo (Figure 4). Additional evidence 
that galaxies moving in the cluster environment 
can produce filaments of cooler gas behind them 
were found by David et al. (1994) or Fabian et al. 
(2001). 

The accurate simulation of dynamic stripping 
in cluster and group elliptical galaxies performed 
by Toniazzo & Schindler (2001), show that the ef- 
fect of dynamic stripping depends on the galaxy 
orbit and the duration of the subsonic and super- 
sonic phases which determine the balance within 
gas losses and accretion onto the galaxy. In gen- 
eral their simulations reveal the formation of com- 
plex and elongated tails behind the galaxy (similar 
to those seen in NGC 4472) and strong decenter- 
ing of X-ray isophotes due to supersonic stripping. 
Such irregular structures are not seen in the NGC 
1399 galactic halo which seems instead rather cir- 
cular {eoalactic = 0.02). 

A very rough estimate of the velocity of NGC 
1399 can be obtained supposing that the displace- 
ment of the optical galaxy, coincident with the X- 
ray centroid, with respect to the galactic halo is 
due to the deceleration of the halo caused by ram 
pressure from the cluster component. In the hy- 
pothesis that the motion is highly subsonic, we 
can treat the different components as homoge- 
neous spheres and calciilatc the initial velocity re- 
quired to produce the observed displacement after 
a time t = 10^ yr due to the ram pressure of the 
cluster halo P ~ PciusterV^ ■ Using the parame- 
ters derived form the bidimensional model (§ 2.5): 
Pgaiaxy/ Pciuster ~ 4, Rgaiaxy ~ 4' and a displace- 
ment AS ii?, we obtain an initial velocity of 

10 km s~^, consistent with the initial assump- 
tion of subsonic motion. Even considering a more 
recent interaction where t = 10* the galaxy motion 
remains highly subsonic (~ 100 km s~^). This re- 
sult is in agreement with the observations, showing 
that the halo morphology is only weakly distorted. 
In conclusion, even though extreme projection ef- 



fects may hide stronger stripping features, ram 
pressure may be causing the galaxy halo displace- 
ment but seems unlikely to be responsible for the 
observed inner halo structure of NGC 1399. 

3.2.4- Cluster Merging effects 

For what concerns the elongated cluster halo 

component, in many cases irregular features in 
the cluster halos seems associated with merger 
events (see for instance Sun et al. 2000; Maz- 
zotta et al. 2000). Simulations of cluster merg- 
ing (Roettiger et al. 1996) show that elongation 
and twisting of the X-ray isophotes are a general 
result of such interactions. In the case of the For- 
nax cluster Drinkwater, Gregg & CoUess (2000) 
found evidence of a subcluster 3 degrees South- 
West of NGC 1399 centered on NGC 1316 (For- 
nax A). They claim that the two components are 
at the early stage of merging; in this case merging 
is unlikely to be the cause of the elongation of the 
cluster X-ray isophotes. To further study the dy- 
namical status of the cluster halo we would need 
temperature maps of the whole cluster, as can be 
obtained with Chandra or XMM. 

It is worth noticing that the X-ray centroid of 
the cluster component is displaced with respect 
to the optical distribution of Fornax cluster galax- 
ies, which tend to concentrate South- West of NGC 
1399 (see Figure 1 from Hilker et al. 1999a and 
Figure 3 from Drinkwater, Gregg & CoUess 2000). 
If the merging picture is confirmed, this may im- 
ply that the galaxies in the two subclusters are 
infalling in the global potential well faster than 
the associated gaseous halos. 

Kim & Fabbiano (1995) compared three X-ray 
dominant galaxies in poor groups: NGC 1399, 
NGC 507 and NGC 5044. They show that even 
if these galaxies have similar properties (tempera- 
tures, cooling core) their X-ray brightness profiles 
are very different. This suggests that the shape of 
the potential and the mass deposition mechanism 
may vary from galaxy to galaxy. If environmental 
effects (tidal interactions, ram pressure stripping) 
are important factors in determining the NGC 
1399 halo structure, they must be considered as 
well in order to explain the observed differences. 
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3.3. The NGC 1404 halo 

In contrast with NGC 1399, NGC 1404 pos- 
sesses a very regular halo, well represented by a 
Beta profile. The X-ray brightness profile falls as 
^ ^-2.08±o.03^ good agreement with the optical 
brightness profile S oc 7-i-9±o i (Forbes et al. 1998). 
This result suggests that the the gas distribution 
of the hot halo of NGC 1404 is produced by the 
same gravitational potential that is binding the 
stellar population, out to r ~ 80" (~ 7 kpc) from 
the galaxy center (Figure 10b). This result is in 
agreement with Matsushita (2001) who found that 
luminosities and temperatures of X-ray compact 
galaxies are well explained by kincmatical heating 
of gas supplied by stellar mass loss. However she 
notice that NGC 1404 is significantly brighter than 
other compact ellipticals and explains this result 
in term of interaction with the ISM. 

Indeed at radii larger than 80" the influence of 
dynamic stripping is evident from the presence of 
an elongated tail on the South-East side of the 
galaxy and the steep surface brightness gradient 
on the opposite side (Figure 3 and 10b). However 
the regularity of the brightness profile and its sim- 
ilarities to the optical profile within 90" (~ 8 kpc) 
implies that either the stripping efficiency is low or 
that it affects significantly only the outer galactic 
regions. It is worth noticing that our brightness 
profile is shallower than expected from simulations 
of dynamically stripped galaxies suggesting that 
the models do not properly describe the influence 
of ram stripping on the inner galactic halo (see 
Toniazzo & Schindler 2001). 

The gas and total mass profiles derived in § 
2.6.2 are in the same range of those of the central 
component of the NGC 1399 halo, within the inner 
80" . This is not surprising since we have seen that 
in NGC 1399 the central component seems to be in 
equilibrium with the stellar population, as it is the 
case for the NGC 1404 halo. This result further 
supports the scenario in which the difference be- 
tween X-ray compact and extended galaxies lies in 
the presence of additional components (i.e. galac- 
tic and cluster components) related to the central 
position of the dominant galaxic^s inside the cluster 
potential well (e.g. Matsushita 2001). 



3.4. X-ray— Radio Interactions 

Signs of interaction between radio jet/lobes 
and the surrounding gaseous medium have been 
reported in powerful radio galaxies (see for in- 
stance Carilli, Perley & Harris 1994 for Cygnus 
A; Bohringer et al. 1993; McNamara et al. 2000; 
Fabian et al. 2000 for evidence of such interaction 
in clusters hosting powerful radio sources). NGC 
1399 is a faint radio galaxy whose radio lobes are 
entirely contained within the optical galaxy. Nev- 
ertheless the structures seen in Figure 20 suggest 
that the same mechanisms at work in more pow- 
erful sources may be responsible for some of the 
galactic halo structures. The residual image in 
Figure 21 shows excess emission on both sides of 
the nucleus and a lack of X-ray emission along 
the radio jets in the inner 15" (^ 1.4 kpc), as if 
the radio jets are displacing the ambient gas. For 
r ~ 25" (~ 2.3 kpc) more X-ray emission is seen 
in the North and South part of the galaxy, aligned 
with the direction of the radio jets. 

Clarke, Harris & Carilli (1997) have studied, by 
means of hydrodynamical simulations, how these 
features are produced during the propagation of 
the radio jet trough the ICM. They show that, for 
a --^ 1 keV gaseous halo and a detector response 
similar to the ROSAT HRI (see their Figure 8), 
we expect to see excess emission located on the 
side of the jets, due to the shocked gas near the 
tail of the radio lobes, as the one seen in Figure 
20 (region A) and in Figure 21. 

In coincidence with the radio lobes, X-ray en- 
hancements or deficits are expected depending on 
the thickness of the shocked gas Ar with respect 
to the lobe radius r. Near the edge of the radio 
lobe, where r > 3.4Ar, the shocked gas density 
is high enough to produce a brightness enhance- 
ment. On the contrary, in the tail of the lobes 
the compressed gas had time to expand, so that 
r < 3.4Ar and the X-ray brightness is diminished, 
producing a 'cavity'. This picture is very similar 
to what we see in our data. Features B and C in 
Figure 20 may suggest respectively the presence 
of a cavity produced by the Southern lobe in the 
ICM and the excess emission associated with the 
edge of the lobe. The alternative hypothesis, that 
the excess is due to the presence of a bow shock 
at the leading end of the lobe, seems unlikely both 
because Carilli and collaborators showed that the 



29 



bow shock is hardly seen at energies < 4 keV and 
because the signatures of bow shocks are usually 
sharper than the one seen in our data, showing 
up as compact bright spots or arcs. The Chan- 
dra image also revealed the presence of a long arc 
South of NGC 1399. The steeper gradient of the 
radio emission on the Western side of the South- 
ern lobe, with respect of the Eastern side, suggests 
that this feature may represent a layer of hot gas 
compressed by the radio emitting plasma. 

Signatures of interaction with the ICM do not 
show up so clearly in the Northern lobe, but this 
could be due to projection effects (see for instance 
Brunetti et al. 2000; Finoguenov & Jones 2000). 
The brightness enhancement near the West side 
of the lobe could be due to the thin layer of com- 
pressed gas surrounding the lobe while the steeper 
radio gradient near region D suggests that the lobe 
is avoiding a region of higher gas density. Indeed 
the crude estimate indicates that the pressure of 
the bubble (§ 2.8) can be up to 4 times higher than 
that of the surrounding ISM. An upper limit of 
k ~ 420 on the ratio of the electron to heavy par- 
ticle energy (or alternatively a lower limit of 0.005 
on the filling factor) can be derived by match- 
ing the radio pressure to the thermal pressure of 
the bubble (Figure 22), in agreement with Killeen, 
Bicknell & Ekers (1988). 

Inverse Compton (IC) scattering of Cosmic 
Background photons by the radio emitting par- 
ticles is often considered as an alternative process 
to explain the X-ray brightness enhancements near 
radio lobes (Harris & Grindlay 1979; Kaneda et 
al. 1995). We have evaluated this possibility for 
the X-ray excess near the edge of the Southern 
radio lobe of NGC 1399, following Feigelson et 
al. (1995), and comparing the minimum energy 
field strength Bme with the magnetic field Bjc 
required to produce a measured ratio of radio to 
X-ray flux. We used the radio flux for the South- 
ern lobe measured by Killeen, Bicknell & Ekers 
(1988), corrected for the fact that just half of the 
lobe contributes to the X-ray excess, and we mea- 
sured the X-ray flux density from the residuals 
shown in Figure 14, assuming that all the excess 
over the 2D model is due to IC scattering. We find 
that it is unlikely that more than a few percent 
of the excess near the edge of the Southern radio 
lobe may be due to IC scattering, even taking into 
account the poor knowledge of the parameters 



involved in the calculations. 

Pressure confinement of the radio lobes (Killeen, 
Bicknell & Ekers 1988) is supported by our data 
(Figure 22). If the inclination of the radio source 
is < 45°, the radio jets remain collimated up to 
the radius within which the high density central 
component dominates (~ 50" , Figure 22) and dif- 
fuse in lobes as soon as they enter in the lower 
density environment. A larger inclination is un- 
likely because the thermal pressure of the cluster 
component (long-dashed line in Figure 22) is of the 
same order of magnitude as the radio lobes pres- 
sure, and would probably not be able to confine 
the radio emitting plasma, if this was extending 
at greater distances. 

Finally, we analized the possibility that the 
power-law component found in the spectral model 
of Allen, di Matteo & Fabian (2000) for NGC 1399 
(see § 3.1) may be due to an active nucleus. Their 
estimated flux (converted to our adopted distance) 
is twice the value expected from our upper limit 
on a possible nuclear source, adopting their best 
fit parameters. Considering that Loewenstein et 
al. (2001) further reduce this upper limit by an 
order of magnitude it appears that the bulk of the 
power-law emission comes from the discrete source 
population revealed by the Chandra data (see Fig- 
ure 25). 

3.5. Discrete sources 

The presence of a population of individual X- 
ray sources in E and SO galaxies was first sug- 
gested by Trinchieri & Fabbiano (1985). Ad- 
ditional evidence was then provided by ROSAT 
and ASCA observations (e.g. Fabbiano, Kim & 
Trinchieri 1984; Kim et al. 1996; Matsumoto et al. 
1997). These X-ray source populations are now 
been detected down to much lower luminosities 
with Chandra (e.g. Sarazin, Irwin & Bregman 
2000, 2001). The HRI data demonstrate the pres- 
ence of such population, possibly of accretion bina- 
ries, also in NGC 1399. This result is confirmed by 
the Chandra data which reveal a large number of 
pointlike sources surrounding the galaxy, many of 
which associated with globular clusters (Angelini, 
Loewenstein & Mushotzky 2001). This population 
of accretion binaries is probably responsible for 
the hard spectral component (Buote 1999; Mat- 
sushita, Ohashi & Makishima 2000; Allen, di Mat- 
teo & Fabian 2000) required to fit the high energy 
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spectrum of NGC 1399. As demonstrated by com- 
paring the 5" resolution Rosat HRl image with the 
10 times sharper Chandra ACIS image, subarcsec- 
ond resolution is essential for the study of X-ray 
source populations in galaxies. 

4. SUMMARY 

We have analized a deep ROSAT HRI obser- 
vation centered on the dominant Fornax clus- 
ter galaxy NGC 1399 and including NGC 1404, 
in conjunction with the archival ROSAT PSPC 
data of the same field. We found an extended 
and asymmetric halo extending on cluster scales 
(r > 90 Kpc). The halo profile is not consistent 
with a simple Beta model but suggests the pres- 
ence of three different components. The combined 
HRI and PSPC data were fitted with a multi- 
component bidimensional model, consisting of: (1) 
a 'central' component (dominating for r < 50") 
centered on the optical galax;y whose distribution 
follows the luminous matter profile; (2) a 'galactic' 
component (50" < r < 400") displaced 1' (5 kpc) 
South- West of NGC 1399 and (3) a 'cluster' ellip- 
soidal component (r > 400") centered ~ 5.6' (31 
kpc) North-East of the galaxy. The HRI image 
also revealed the existence of filamentary struc- 
tures and cavities in the galactic halo, due to den- 
sity fluctuations in the hot gas. The presence of 
these features is confirmed by a preliminary analy- 
sis of Chandra data. The large scale surface bright- 
ness distribution can be explained if the galaxy 
hosts a large dark halo with different dark matter 
distributions on galactic and cluster scales (Saglia 
et al. 2000; Ikebe et al. 1996). Alternative models 
that explain the GCs abundance and optical ve- 
locity dispersion profiles through tidal interactions 
with Fornax cluster galaxies (Kissler-Patig et al. 
1999; Napolitano, Arnaboldi & Capaccioli 2001, in 
preparation) are compatible with our data only if 
such encounters produce a significative flattening 
of the outer gas distribution. Tidal interactions 
may also explain the presence of the density fluc- 
tuations in the galactic halo (D'Ercole, Recchi & 
Ciotti 2000). Ram pressure from the cluster halo 
(e.g. Irwin & Sarazin 1996; David et al. 1994; 
Fabian et al. 2001) is able to account for the de- 
centering of the optical galaxy with respect to the 
galactic halo if NGC 1399 is moving subsonically 
(10-100 km s~^) in the cluster potential. The dis- 
placement of the center of the cluster component 



with respect to the nearest Fornax galaxies may 
suggest that the cluster is not relaxed and may be 
undergoing a merger event. 

We do not detect the X-ray counterpart of the 
nuclear radio source (Killeen, Bicknell & Ekers 
1988) and pose an upper limit of = 3.9 x 10^^ 
erg s"-^ on its luminosity in the 0.1-2.4 keV band. 
In agreement with Killeen and collaborators, our 
data support thermal confinement of the radio 
emitting plasma. We found X-ray emission en- 
hancements aligned with the radio jets that may 
be due to shocked gas and X-ray 'holes' and en- 
hancements coincident with the position of the ra- 
dio lobes. The latter are consistent with a sce- 
nario in which the hot gas is displaced by the ra- 
dio plasma pressure (e.g. Clarke, Harris & Carilli 
1997), while alternative models of Inverse Comp- 
ton scattering of CMB photons fail to account for 
the observed excesses. 

NGC 1404 possess a very symmetric halo within 
8 kpc well fitted by a simple Beta model. The 
agreement with the optical profile and the similar- 
ity with the 'central' component of the NGC 1399 
halo suggest that the ISM distribution of NGC 
1404 is produced by the same gravitational po- 
tential that is binding the stellar population. At 
larger radii the influence of dynamic stripping is 
evident from the presence of an elongated tail on 
the South-East side of the galaxy and the steep 
surface brightness gradient on the opposite side. 

We detected 43 discrete sources in our field. 
Their flux distribution is consistent at the 2(7 level 
with the number expected from X-ray background 
counts. However, there is a significant excess 
'peak' on the central galaxy in the spatial distri- 
bution of these sources, suggesting a population of 
galaxian X-ray sources in NGC 1399. This is con- 
firmed by the Chandra data, that reveal a large 
number of point-like sources associated with the 
optical galaxy. This sources are likely to explain 
the hard component found by ROSAT and ASCA 
spectral analysis (Buote 1999; Matsushita, Ohashi 
& Makishima 2000; Allen, di Matteo & Fabian 
2000). 
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